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Abstract. We present the results of the first extensive mid-infrared (IR) imaging survey of the p Ophiuchi 
embedded cluster, performed with the ISOCAM camera on board the ISO satellite. The main p Ophiuchi molecular 
cloud L1688, as well as the two secondary clouds L1689N and L1689S, have been completely surveyed for point 
sources at 6.7 nm and 14.3 pm. A total of 425 sources are detected in ~ 0.7 deg'^, including 16 Class I, 123 
Glass II, and 77 Class III young stellar objects (YSOs). Essentially all of the mid-IR sources coincide with near-IR 
sources, but a large proportion of them are recognized for the first time as YSOs. Our dual-wavelength survey 
allows us to identify essentially all the YSOs with IR excess in the embedded cluster down to Fi, ~ 10-15 mjy. 
It more than doubles the known population of Glass II YSOs and represents the most complete census to date 
of newly formed stars in the p Ophiuchi central region. There are, however, reasons to believe that several tens 
of Class III YSOs remain to be identified below L* ~ 0.2 Lq. The mid-IR luminosities of most (~ 65%) Class II 
objects are consistent with emission from purely passive circumstellar disks. The stellar luminosity function of 
the complete sample of Glass II YSOs is derived with a good accuracy down to Li, ~ 0.03 L0. It is basically fiat 
(in logarithmic units) below L* ~ 2 Lq, exhibits a possible local maximum at L* ~ 1.5 Lq, and sharply falls 
off at higher luminosities. A modeling of the luminosity function, using available pre-main sequence tracks and 
plausible star formation histories, allows us to derive the mass distribution of the Class II YSOs which arguably 
reflects the initial mass function (IMF) of the embedded cluster. After correction for the presence of unresolved 
binary systems, we estimate that the IMF in p Ophiuchi is well described by a two-component power law with a 
low-mass index of —0.35 ± 0.25, a high-mass index of —1.7 (to be compared with the Salpeter value of —1.35), 
and a break occurring at Mflat ~ 0.55 ± 0.25 Mq. This IMF is flat with no evidence for a low-mass cutoff down 
to at least ~ 0.06 Mq. 

Key words. Stars: formation - Stars: low-mass, brown dwarfs - Stars: luminosity function, mass function - Stars: 
pre-main sequence - ISM: individual objects: p Ophiuchi cloud 
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1. Introduction 

Recent observations suggest that most stars in our Galaxy 
and other galaxies form in compact clusters. In particular, 
near-IR imaging surveys of nearby molecular cloud com- 
plexes have shown that the star formation activity is typ- 
ically concentrated within a few rich clusters associated 
with massive dense cores which constitute only a small 
fraction of the total gas mass available (e.g. Lada, Strom, 
& Myers |1993| , Zinnecker, McCaughrean, & Wilking |l99|) . 
These embedded clusters comprise various types of young 
stellar objects (YSOs) - from still collapsing protostars to 
young main sequence stars - but are usually dominated 
in number by low-mass pre-main sequence (PMS) stars, 
i.e., T Tauri stars. Young clusters provide excellent labo- 
ratories for studying the formation and early evolution of 
stars through the observational analysis of large, geneti- 
cally homogeneous samples of embedded YSOs. Two key 
characteristics of these young stellar populations are their 
luminosity distribution and their mass spectrum, which 
give important observational constraints on the stellar ini- 
tial mass function (IMF) in the Galaxy. Observations of 
young embedded clusters can also help us understand pos- 
sible links between parent cloud properties and the result- 
ing stellar masses. They are however hampered by: (1) 
dust extinction from the parent cloud which hides most of 
the newly formed stars at optical wavelengths, (2) the dif- 
ficulty to recognize the nature of individual sources (e.g. 
protostars, T Tauri stars, or background sources), (3) the 
youth (and thus poorly known intrinsic properties) of most 
cluster members. 

Owing to these difficulties, the census of embedded 
YSOs provided by IRAS and near-IR studies is far from 
complete even in the nearest clouds (e.g. Wilking, Lada, 
& Young |1989| - hereafter WLY89). Thanks to its high 



sensitivity and good spatial resolution in the mid-I R, the 
ISOCAM ca mera o n board ISO (Cesarsky et al 



1996 



Kessler et al. |1996D was a powerful tool to achieve more 
complete surveys for YSOs in all major n earby star- 
forming regions 
al 



1999; Persi et al. 200C) 



see Nordh et al. 1996. 1998; Olofsson et 



The nearby p Ophiuchi cloud is one of the most ac- 
tively studied sites of low-mass star formation. Its cen- 
tral region harbors a rich embedded cluster with about 
100 members recognized prior to the present work (e.g. 
WLY89, Casanova et al. 1995). While a dispersed pop- 
ulation of optically visible young stars, associated with 
the Upper-Scorpius OB association, has a typical age of 
several million years (Myr) (e.g. Preibisch & Zinnecker 



1999), the central embedded cluster is recognized as one 



of the youngest clusters known with an estimated age 
on the order of ~ 0.3-1 Myr (e.g. WLY89, Greene & 



Meyer 1995, Luhman & Rieke 1999). This young cluster 
has been extensively studied at wavelengths ranging from 
the X-ray to the radio band. The satellites Einstein and 
ROSAT have revealed ^ 70 highly variable X-ray sources 
associated with magnetically-active young stars, includ- 
ing deeply embedded protostellar sources (Montmerle et 



al. 1983, Casanova et al. 1995, Grosso et al. 2000). In 
the near-IR, the cloud has been deeply surveyed from the 
ground using large-format arrays (Greene & Young 1992, 
Comeron et al. 1993, Strom et al. 1995, Barsony et al. 



1997). Unfortunately, due to difficulties in discriminating 



between background sources and embedded YSOs with- 
out performing time-consuming mid-IR photometry (e.g. 
Greene et al. 1994) or near-IR spectroscopy (e.g. Greene & 
Lada 199(:; Luhman & Rieke 1999| ), these recent near-IR 
surveys have only partially increased the number of clas- 
sified, recognized members. Finally, while only relatively 
poor angular resolution IRAS data are available so far in 
the far-IR (e.g. WLY89), deep imaging surveys at an angu- 
lar resolution of 10 — 15" or better exist a t (sub )millimeter 



wavelengths (e.g., Motte , And re, & Neri 1998| - hereafter 



MAN98 
It 



-; Wilson et al 
thanks to the 



199C) 



IS thanks to ttic illuminating example of the 
p Ophiuchi embedded cluster that the now widely used 
empirical classification of YSOs was originally introduced. 
Three IR classes were initially distinguished based on 
the shapes of the observed spectral energy distributions 
(SEDs) between ^ 2 /im and ~ 25-100 p,\a (Lada & 
Wilking |1984 WLY89). Objects with rising SEDs in this 
wavelength range were classified as Class I, sources with 
SEDs broader than blackbodies but decreasing longward 
of ~ 2 /im as Class II, and sources with SEDs consis- 
tent with (or only slightly broader than) reddened stel- 
lar blackbodies as Class III. These morphological SED 
classes are interpreted in terms of an evolutionary se- 
quence from (evolved) protostars (Class I), to T Tauri 
stars with optically thick IR circumstellar disks (Class II), 
to weak T Tauri stars with at most opticall y thi n disks 
(Class III) (La da |l987t Adams, Lada, & Shu jl987t Andre 
& Montmerle |l994| - hereafter AM94). A fourth class 
(Class 0) was subsequen tly in troduced by Andre, Ward- 
Thompson, & Barsony ( 1993 ) to accommodate the dis- 
covery in the radio range of cold sources with large sub- 
millimeter to bolometric luminosity ratios and powerful 
jet-like outfl ows, such as VLA 1 623 i n p Oph A (e.g. 
Andre et al. 1990; Bontemps et al. 1996| ). Class objects, 
which have measured circumstellar envelope masses larger 
than their inferred central stellar masses, are interpreted 
as young protostars at the beginning of the main accretion 
phase (e.g. Andre, Ward-Thompson, & Barsony 2000). 
The fact that the p Ophiuchi central region contains at 
least two Class protostars (Andre, Ward-Thompson, & 
Barsony 1993| )), as well as numerous ( ^ 60) pre- stellar 
condensations (MAN98), demonstrates that it is still ac- 
tively forming stars at the present time. 

The distance to the p Ophiuchi cloud is somewhat un- 
certain. Usually, a value of 160 pc is adopted (e.g. Chini 
1981). However, recent Hipparcos results on the Upper- 



Scorpius OB association (de Zeeuw et al. 1999) provide a 
reasonably accurate estimate of 145 ± 2 pc for the average 
distance to the stars of the OB association. The embed- 
ded cluster is located at the inner edge of the molecular 
comp lex on the outskirts of the OB association (e.g. de 
Geus 198£), and not very far, in projection, from the as- 
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sociation center (less than 4° apart which corresponds to 
~ 10 pc; see Fig. 9 of de Zeeuw et al. 199!;). In this pa- 
per, we therefore adopt a distance of d = 140 ± 10 pc for 
the p Ophiuchi IR cluster which corresponds to a distance 
modulus 5xlogio(rf/10pc) = 5.73 ± 0.15. 

The layout of the paper is as follows. Sect. 2 gives ob- 
servational details (§ 2.2) and describes the way the data 
have been reduced to obtain mid-IR images and extract 
point-sources (§ 2.3) along with the photometric uncer- 
tainties and the sensitivity levels (§ 2.4). In Sect. 3, the 
identification of detected sources is discussed (§ 3.1) and 
the selection of a new population of 123 Class II YSOs, 
as well as 16 Class I and 77 Class III YSOs, is described 
(§ 3.2-3.5). We then derive luminosities for these YSOs and 
build the corresponding luminosity functions in Sect. 4. In 
Sect. 5, we model the luminosity functions for Class II and 
Class III YSOs in terms of the underlying mass function 
and star formation history. In Sect. 6, we discuss the re- 
sulting constraints on the IMF of the p Ophiuchi cloud 
down to ^ 0.06 Mq, (§ 6.1-6.2) as well as related implica- 
tions (§ 6.3-6.5). 



2. Observations and data reduction 

2.1. Region surveyed by ISOCAM 

The survey encompasses the p Ophiuchi central region as- 
sociated with the prominent dark cloud L1688, as well as 
the two subsidiary sites L1689N and L1689S. The L1688 
field is a ~ 45' x 45' square, while the L1689N and L1689S 
fields each cover an area of 16.5' x 16' (see Fig. 1). Most 
previously known members of the p Ophiuchi cluster lie 
within these fields. In particular, this is the case for 94 of 
a total of 113 recognized members from WLY89, AM94, 



Greene et al. (1994), and MAN98. The known young stars 
which lie outside the boundaries of our survey are mostly 
optically visible, weak-line or post T Tauri stars (belong- 
ing to Class III) spread over a large area on the o utskirts 
of the molecular complex (e.g. Martm et al. 199S). 



2.2. Observational details 

The mapping was performed in the raster mode of 
ISOCAM in which the mid-IR 32 x 32 pixel array imaged 
the sky at consecutive positions along a series of scans par- 
allel to the right-ascension axis. The offset between consec- 
utive array positions along each scan (Aa) was 15 pixels, 
while the offset between two scans (A(5) was 26 pixels. 
Each set of scans was then co-added and combined into 
a single raster image. The final image of the L1688 field 
(Fig. 1 and Abergel et al. 1996) actually results from the 
combination of six separate rasters. A pixel field of view 
of 6" was used for four of these rasters, but smaller 3" 
pixels were employed for the other two rasters in order to 
avoid saturating the array on the brightest sources of the 
cluster. The L1689N and L1689S fields were imaged with 
one raster each using 6" pixels. 



In order to avoid saturation, the individual readout 
time for the L1688 rasters was set to tint = 0.28 s. About 
55 of these readouts (i.e. an integration time of ^ 15 s) 
were performed per sky position. Thanks to the half-frame 
overlap between subsequent individual images, each sky 
position was observed twice, yielding an effective total in- 
tegration time of - 30 s. For L1689N and L1689S it was 
possible to use ijnt = 2.1s, and about 15 readouts were 
performed per sky position with the same half-frame over- 
lap, giving an integration time per sky position of ~ 60 s. 
A total of 1104 individual images were necessary to mo- 
saic the L1688 field, and an additional 60 images each were 
used to map L1689N and L1689S. 

All three fields were mapped in two broad-band filters 
of ISOCAM: LW2 (5-8.5/im) and LW3 (12-18/ini). These 
filters are approximately centered on two minima of the 
interstellar extinction curve and are situated apart from 
the silicate absorption bands (at roughly 10 and 18 p.m). 
However, they include most of the Unidentified Infrared 
Bands (UIBs, likely due to PAH-like molecules) which con- 
stitute a major source of backgroun d emis sion toward star- 
forming clouds (e.g. Bernard et al. 1993, Boulanger et al. 



1996 ). The ISOCAM central wavelengths adopted here for 



LW2 and LW3 are 6.7 /im and 14.3 /im respectively. 

2.3. Image processing, source extraction, and 
photometry 

Each raster consists of a temporal series of individual 
integration frames (i.e. of 32 x 32 pixel images) which 
was reduced using the CAM Interactive Analysis software 
(CIA)Q. We have subtracted the best dark current from 
the ISOCAM calibration library, and as a second step we 
improved it with a second order correction using a EFT 
thresholding method (Starck et al. 1999). Cosmic-ray hits 
were detected and masked using the mu lti-res olution me- 
dian transform algorithm (Starck et al. 1996 ). The tran- 
sients in the time history of each pixel due to detector 
memory effects were correct ed wi th the inversion method 
described in Abergel et al. ( 1996| ). The images were then 
flat-fielded with a fiat image obtained from the observa- 
tions themselves. Since these various corrections applied to 
the images are not perfect, the extraction of faint sources 
from the images is a difficult task. We have developed 
an interactive IDL point-source detection and photometry 
program for raster observations which works in the CIA 
environment. This program helps to discriminate between 
astronomical sources and remaining low-level glitches or 
ghosts due to strong transients (see also Nordh et al. 



1996; Kaas et al. 2001). The fiuxes of the detected sources 



were estimated from the series of flux measurements made 
in the individual images (usually 2 to 4 individual im- 
ages cover each source) which were obtained from classical 



^ CIA is a joint development by the ESA Astrophysics 
Division and the ISOCAM consortium led by the ISOCAM 
PI, C. Cesarsky, Direction des Sciences de la Matiere, C.E.A., 
France. 
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Fig. 1. Sky map of the p Ophiuchi clouds L1688, L1689N, and L1689S showing the spatial distribution of the ISOCAM 
sources. The filled symbols correspond to "red" YSOs (i.e. Class I and Class II YSOs, see Sect. 3), while the open 
circles mark the "blue" sources. The surface areas of these symbols are proportional to the 6.7/zm flux densities. The 
open stars mark four young early-type stars (SR3, SI, WL22, WL16 from right to left; see Sect. 3.1). The contours 
show the L1688 molecular cloud as mapped in CS(2-1) by Liseau et al. (1995). The triangles refer to the peaks of 
the DCO^ dense cores of Loren et al. (1990). The figure emphasizes the presence of four sub-clusters Oph A, B, EF 
(between Oph E and Oph F), and L1689S (named after the associated dense cores). 



aperture photometry. The emission was integrated in a sky 
aperture, the background emission subtracted, and finally 
an appropriate aperture correction was applied based on 
observed point-spread functions available in the ISOCAM 
calibration library. In practice, the radius of the aperture 
used was 9" (i.e., 3 and 1.5 pixels for a pixel size of 3" 
and 6", respectively). For the weakest sources, however, 
we reduced the aperture radius to 4.5" (i.e. 1.5 pixels for 
a pixel size of 3"), in order to improve the signal-to- noise 
ratio. Finally, we applied the following conversion factors: 
2.33 and 1.97 ADU/gain/s/mJy for LW2 and LW3 re- 
spectively (from in-orbit latest calibration - Blommacrt 
199S ) . These calibration factors are strictly valid only for 
sources with a flat SED (F^ ~ Here, a small but 

significant ( <L 1 %) color correction needs to be applied 
to the bluest sources, recognized as Class III YSOs in 
Sect. 3 below. For these sources, the conversion factors 
quoted above were divided by 1.05 for LW2 and 1.02 for 
LW3 to account for the color effect. The 212 ISOCAM 
sources recognized as cluster members (see Sect. 3 below) 
are listed in Table 1 (available only in electronic form 
at http://cdsweb.u-strasbg.fr/) with their J2000 coordi- 



ties (see Sect. 2.4), as well as the corresponding near-IR 
identifications. 



2.4. Photometric uncertainties and point-source 
sensitivity 

The uncertainties on the final photometric measurements 
result from systematic errors due to uncertainties in the 
absolute calibration and the aperture correction factors, 
and from random errors associated with the fiat-fielding 
noise, the statistical noise in the raw data, the noise due 
to remaining low-level glitches, and the imperfect correc- 
tion for the transient behavior of the detectors. The in- 
orbit absolute calibration has been verified to be correct 
to within 5 % (Blommaert 1998), and we estimate that 



nates, their fiux densities and associated rms uncertain- 



the maximum systematic error on the aperture correction 
is ~ 10 % (by comparing theoretical and observed point- 
spread functions). The maximum systematic error on our 
photometry is thus ^ 15 %. The magnitudes of the random 
errors were directly estimated from the data by measuring 
both a "temporal" noise (noise in the temporal sequence of 
individual integrations) and a "spatial" noise (due to im- 
perfect flat-fielding and/or spatial structures in the local 
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mid-IR background emission) for each source in the auto- 
matic detection procedure. The temporal noise was com- 
puted as the standard deviation of the individual aperture 
measurements divided by the square root of the number 
of measurements. The spatial noise was estimated as the 
standard deviation around the mean background (linear 
combination of the median and the mean of the pixels op- 
timized for the source flux estimates) in the immediate 
vicinity of each source. 
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Fig. 2. (a) Distribution of the F^ '^ flux densities for the 
425 sources detected by ISOCAM. The dark histogram 
shows the flux distribution of the cluster members (i.e., 
YSOs) discussed in Sect. 3, while the light histogram cor- 
responds to the remaining sources which are likely dom- 
inated by background Galactic objects, (b) Same as (a) 
for the F^"^-^ flux densities. The dark histogram shows 
the flux distribution for the "red" ISOCAM sources (i.e., 
Class I and II YSOs), while the light histogram comprises 
only the "blue" ISOCAM sources not associated with rec- 
ognized YSOs (see Sect. 3 and Fig. 3). The adopted com- 
pleteness levels equal to 10 and 15mJy for F^ '' and F^'^-^ 
respectively are shown as vertical dotted lines. The solid 
and dashed curves show the expected number of Galactic 
sources in each bin according to the model of Wainscoat 
et al. (1992) for Ay = and Ay = 10 respectively. 



The sensitivity limit of the survey was estimated by 
calculating the average value of the quadratic sum of the 
temporal and spatial noises measured on the weakest de- 
tected sources. The total rms flux uncertainty found in 
this way, atot = {<^tenip + ^^Ipatf^^ is ^6.7 = 2.2 mJy at 
6.7 /im and cri4.3 — 4.1 mJy at 14.3 /im, 75 % of which 
is due to the spatial noise component. The large contri- 
bution of the spatial noise originates in the highly struc- 
tured diffuse mid-IR emission from the ambient molecular 
cloud itself (see Abergel et al. 



1996). Figure 2 displays 



the distributions of fluxes at 6.7 /xm and 14.3 ixm for all 
the detected ISOCAM sources. We used the Wainscoat 
et al. ( 1992] ) Galactic model of the mid-IR point source 
sky to estimate the expected number of foreground and 
background sources up to a distance of 20 kpc. The model 
predictions are shown by solid and dashed curves in Fig. 2 
for cloud extinctions of Ay — and Ay — 10, respec- 



seen that the flux histograms of the ISOCAM sources not 
associated with YSOs (light shading in Fig. 2) are remark- 
ably similar in shape to the model distributions down to 
^ 6 mJy at 6.7 //m and ^ 10 mJy at 14.3 /im. These flux 
densities can be used to estimate the completeness level 
of our observations which is not uniform over the spa- 
tial extent of the survey. The histograms with light (grey) 
shading in Fig. 2 are dominated by background sources 
preferentially located in low-noise regions (i.e., outside the 
crowded central part of L1688), where the total rms flux 
uncertainty is ~ 2.0 mJy at 6.7 /im and 3.5 mJy at 
14.3 /tm. The effective completeness level in these regions 
is thus ~ icFtot, where atot is the totalQ flux uncertainty 
(see above). However, most of the YSOs are located in re- 
gions where the noise is somewhat larger. The largest rms 
noise is reached in the Oph A core area (see Fig. 1), where 
o'tot ~ 3.4 mJy at 6.7 /tm and atot ~ 5.0 mJy at 14.3 /im. 
Therefore, we conservatively estimate the completeness 
levels of the global ISOCAM survey to be ~ 10 mJy at 
6.7 /tm and 15 mJy at 14.3 /im. 

Finally, we note that the Ay = 10 model curve in 
Fig. 2b accounts for essentially all the "blue" sources de- 
tected at 14.3 fjim and not associated with known YSOs. 
At 6.7 /im, the predictions of the Wainscoat et al. model 
suggest that there might still be a slight excess of ^ 30 
unidentifled sources belonging to the cloud (Fig. 2a). 



3. Identification and nature of the mid-IR sources 

3.1. Statistics of detections 

Within the ^ 0.7 square degree imaged by ISOCAM, a to- 
tal of 425 sources has been identifled, among which 211 are 
seen at both 6.7 /tm and 14.3 /im. The spatial distribution 
of the sources is shown in Fig. 1, where the "red" sources 
[those with logio{F^^-^ / F^^-'^) > -0.2, see below] are in- 
dicated as filled circles. These "red" sources appear to be 
clustered into four main groupings: three sub-clusters in 
L1688, i.e., Oph A (West), Oph B (North- East) Oph EF 
(South) (see also Strom, Kepner, & Strom 1995), as well 
as a new sub-cluster in L1689S. 

Four bright embedded stars (SI, SR3, WL16, WL22) 
are spatially resolved by ISOCAM in both filters. Their 
extended mid-IR emission is most likely due to PAH- 
like molecules excited by relatively strong far-ultraviolet 
(FUV) radiation fields^. These bright sources are dis- 
played as open star symbols in Fig. 1 (SR3, SI, WL22, 
WL16 from right to left). 



tively (see Kaas et al. P00l| for more details). It can be 



^ The detection level is more directly related to the temporal 
noise than to the spatial noise. In terms of atemp, the effective 
completeness level is ~ 12atemp. 

^ SI and SR3 are optically- visible stars of spectral types B3 
and B7, respectively, which clearly emit enough UV photons to 
excite PAHs. The other two objects, WL16 and WL22 , show 
PAH features in their mid-IR spectra (e.g. Moore et al. 1998 ), 
and may be young embedded early-type (i.e., early A or late 
B) stars. 
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A total of 89 previously classified YSOs lie within the 
area of the present survey: 2 Class 0, 69 Class I/II, and 18 
Class III YSOs. ISOCAM detected 84 of these 89 YSOs 
(94 %): 97 % of the Class I/IIs, 94 % of the Class Ills, 
and none of the Class Os. The two undetected Class I/IIs 
(GY256, GY257), and the undetected Class III (IRS50), 
are located very close (~ 10"-15") to bright sources (WL6 
and IRS48, respectively), which may account for their 
non-detection. The two Class objects (VLA 1623 and 
IRAS 16293—2422) are deeply embedded within massive, 
cold circumstellar envelopes which are probably opaque 
at 6.7 /i m an d 14.3 /im (see Andre, Ward- Thompson, & 
Barsony 1993| )) and too weak to be detected by ISOCAM. 



3.2. A population of new YSOs with mid-1 R excess 

As can be seen in Fig. 3, the mid-IR regime is ideal to 
detect and characterize the excess emission due to circum- 
stellar disks around young stars. In Fig. 3a, the 207 sources 
detected at 6.7 ^m and 14.3 ^m (excluding SI, SR3, 
WL16, WL22) are shown in a diagram which displays the 
logarithmic flux ratio logio(i^J*'^/i^^''') against the mid- 
IR flux Fl'^-^ . On the right axis of the diagram, the flux ra- 
tio F^'^-^ / F^-"^ has been converted into a classical IR spec- 
tral index, aiR = dlogio(AFA)/dlogio(A) (e.g. WLY89), 
calculated between 6.7 /im and 14.3 fxm, i.e., a[^^'^. In 
this diagram, two groups of sources can clearly be distin- 
guished. The lower ratio group has logio(i^J^ '^/i^^' '') ^ 
—0.7, i.e., a spectral index on the order of Qfjj^^* ~ —3.0, 
which is the value expected for simple photospheric black- 
body emission in the Rayleigh- Jeans regime. The disper- 
sion around Qfjj^'^'* ~ —3.0 is obviously larger for weaker 
sources. This is mainly due to increasing photometric un- 
certainty with decreasing flux. The higher ratio group con- 
sists of "red" sourc es de fined by \ogw{F^^-^ / F^ '^) > -0.2 



(cf. Nordh et al. |l996D 
-0.2 < logio(Fj4-3/F6-7) 



Most of these actually have 
< 0.3, i.e., -1.6 



0, 



which is typical of classical T Tauri stars (hereafter 
CTTS). This range of ajj^^^ roughly delineates the do- 
main of Class II YSOs (e.g. Adams et al. 1987, AM94, 
Greene et al. 1994), and is usually interpreted in terms 
of optical ly thi ck circumstellar disk emission (e.g. Lada 



& Adams [19921) . In Fig. 3, the new ISOCAM sources are 
distinguished from the previously known cluster members 
by different symbols. One can see in Fig. 3a that all the 
previously known Class II sources but one lie above the di- 
viding line for red sources, while all the previously known 
Class III sources but two lie below it. 

The object WL19, classified as a Class I YSO by 
WLY89 and as a reddened Class II by AM94 (see also 
Lada & Wilking 1984), is here found to be a "blue" 
source in the mid-IR range (a[j^^^ = —1.8 ± 0.2). This 
object may correspond to a luminous Class III s tar lo - 
cated behind the cloud (see also Comeron et al. 1993|) . 
Although GY12 formally has a "red" mid-IR spectral in- 
dex (a^R = —0.5 ± 0.5) here, we still consider it as 



color is highly uncertain since GY12 is only marginaly re- 
solved at 14.3 /im from its bright Class I neighbor GSS30.) 
Finally the Class III source DoAr21 has a borderline mid- 
IR spectral index {al^'^'^ — —1.5 ±0.1) but is kept as a 
Class III object since its color between 2.2 /im and 14.3 /im 
corresponds to ajj^^"' = —2.0 ± 0.1. 

A total of 71 sources are identified for the first time as 
mid-IR excess objects in Fig. 3a. These new "red" sources 
are most likely all embedded YSOs, i.e., members of the 
p Ophiuchi cluster. Since dust extinction is roughly the 
same at 6.7 /im and 14.3 /im (Rieke & lebofsky 1985; Lutz 
1999 ), any background source should be intrinsically red 
in order to contaminate the sample of red Y SOs. B ased 
on the Galactic model by Wainscoat et al. ( 1992 ), the 
vast majority of background objects should appear "blue" 
(cf. Fig. 2). Seven background giant stars (GY45, GY65, 
VSSG6, GY232, GY351, GY411, and GY453 - cf. Luhman 
& Rieke 1999D, a nd a known foreground dwarf (HD148352 
- Garrison 1967) are detected, which all have blue mid-IR 



colors (ajj^ 



-2.7, -2.0, -3.1, -3.3, -2.9, -3.2, -3.1, 



and —3.1, respectively). 

The newly identified cloud members nicely extend the 
previously known Class I/II population toward low IR 
fluxes. While previous studies could only identify sources 
with F^'^-^ ^ 100 mJy, the present census is complete for 
objects down to F^'^^^ ~ 15 mJy. Altogether, a sample of 
139 Class I/Class II YSOs is identified, of which 71 are 
new members. The present survey has thus allowed us to 
more than double the number of recognized YSOs with 
circumstellar IR excess in the p Ophiuchi cloud. 

The 139 red YSOs are listed in Table 2 (Class I YSOs) 
and Table 3 (Class II YSOs) by decreasing order of cti^^^- 
Twelve of them (last entries of Table 3) are completely new 
sources with respect to published IR surveys. 

3.3. Class I versus Class II YSOs 

In Fig. 3, the YSOs classi fied a s Class I and Class II by 
AM94 and Greene et al. (1994) are shown as filled stars 
and filled circles, respectively. In the \ogio{F^'^-^ / F^ "^) 
versus F^'^-^ diagram, there is no clear color gap between 
the two classes of objects, even though the Class I YSOs 
tend to lie in the upper part of the "red" group (see 
Fig. 3a). Since extinction has a negligible effect on the 
mid-IR ratio Fj'* '^/i^^-^, this suggests that Class I and 
Class II objects have fairly similar intrinsic colors between 
6.7 /im and 14.3 /im. 

The classical IR spectral index calculated from 2 /im to 
10 /im (or 25 pm) (e.g. Lada & Wilking |l98| and WLY89) 
appears to provide a better way of discriminating between 
envelope-dominated Class I YSOs and disk-dominated 
Class II sources (see Fig. 3b). In particular, millimeter 
continuum mapping of optically thin circumstellar dust 
emission confirms that, apart from a few important ex- 
ceptions (e.g., WL22, WL16, WL17, IRS37, IRS47), the 
p Ophiuchi objects selected on the basis of Qfjj^^° <L are 



a Class ill object (cf. Greene et al. 1994| ). (The mid-IR indeed Class I protostars surrounded by spheroidal en 
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Fig. 3. (a) Logarithmic ratio of the ISOCAM fluxes, logwiF^^ '^ / F^-"^), versus F^^-^ displayed for the 207 ISOCAM 
sources seen in both filters (excluding the 4 early-type young stars). The purely ISOCAM color Fj^ '^/F^-^ is converted 
into an IR spectral index, a[^^'^ , on the right axis, and can be used to identify sources with mid-IR excesses (above the 
dashed line). The completeness limit derived in Sect. 2.4 is displayed as a dotted line, (b) Mid-IR color-color diagram 
expressed as a spectral index diagram ajjT^ versus ajjT^'* for 175 ISOCAM sources also detected in the K band (mg.y 



and mi4.3 on the right and upper axes are calculated as in Persi et al. 2000). The inclined heavy dashed line marks the 



limit to the right of which there is significant mid-IR excess as in (a). The two vertical dashed lines mark the formal 
boundary between Class II, flat-spectrum, and Class I sources (see text for details), (c) Classical near-IR color-color 
diagram, J-H versus H-K, for the 152 ISOCAM sources with JHK photometry. 



velopes (AM94, MAN98). 

As expected, the previously known Class I YSOs are well 
concentrated in the upper-right part of the aj^^ 

to the ajpj^ 
and Greene et al 
in this diagram at 



diagram of Fig. 3b (where ^ 

index used in previous studies 



versus 

and ajj^^ are close 



1994) 

T^- 



e.g., AM94 
There is also a hint of two gaps 
ctjj^^^ ^ 0.5 and ajj^^^ ^ 0.0, which 
roughly bracket the regime of flat-spectrum sources as de- 
fined by Greene et al. (1994). These may represent a dis- 
tinct population of transition objects between Class I and 
Class II (e.g. Calvet et al. |1994D . 

with ctTD^"' > 0.55 as 

2-14 
*IR 



Here, we thus consider sources with ajj^'^'* > 
Class I YSOs (Table 2), sources with -0.05 < ajj^"' < 
0.55 as candidate flat-spectrum objects (see Table 3), and 
sources with a^-i4 <• 9.55 as Class II YSOs (Table 3). 
These limiting indices are displayed in Fig. 3b. In the fol- 
lowing, the candidate flat-spectrum sources will be treated 
as Class II YSOs. 



3.4. Mid-IR excess versus near-IR excess 

Most of the ISOCAM sources (e.g. 90 % of the Class II 
sources) were also detected in the near-IR (JHK) survey 
of Barsony et al. ( 1997 ). Comparison between Fig. 3b and 
Fig. 3c illustrates the advantage of mid-IR measurements 
for selecting sources with intrinsic circumstellar IR ex- 
cesses. While the red and blue groups of Fig. 3a are well 
separated in the mid-IR diagram of Fig. 3b, they blend 
together in the near-IR diagram of Fig. 3c. 

Figure 3c also shows that most of the ISOCAM- 
selected YSOs lie within the reddening band associated 
wit h the intrinsic locus of CTTSs as derived by Meyer et 
al. ( 1997 ) in Taurus. The few exceptions, which lie to the 
right of the reddening band, correspond to Class I and 
flat-spectrum YSOs. 
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3.5. The Class III YSO census in p Ophiuchi 

While Class I/II YSOs are easily recognized in the mid- 
IR range thanks to their strong IR excesses, Class III ob- 
jects are difficult to identify without deep X-ray and/or 
radio centimeter continuum observations. We have used 
the ROSAT X-ray surveys of Casanova et al. ( |1995D and 
Grosso et al. ( 2000] ), along with the VLA radio surveys 
by, e.g. , And re, Montmerle, & Feigclson ( |1987D and Stine 
et al. ( 1988 ) to build up a sample of bona- fide Class III 
YSOs covered by the present survey. With the additional 
Class III candidate WL19 (see Sect. 3.2 above), there are 
38 such YSOs which are hsted in order of decreasing a^^^^ 
in Table 4. 

This Class III sample is unfortunately not as com- 
plete as the Class I and Class II samples discussed above. 
According to Grosso et al. ( ^000 ), the number of Class Ills 
may be roughly as large as the number of Class lis: above 
their typical X-ray detection limit of Lx ~ 3 x 10^^ erg s~^ 
(corresponding to L^, ^ 0.3 Lq - see Fig. 7 of Grosso et 
al.), they found a Class Ill/Class II number ratio of 19/22 
in the R.OSAT-HRI/ISO-ISOCAM overlapping survey 
area. If this ratio is representative of the complete popu- 
lation of young stars in p Ophiuchi, the total number of 
Class lis found here (123 objects) suggests that as many 
as 106 Class Ills may be present in the cluster down to 

^ 0.03 Lq (our completeness level for Class lis, see 
Sect. 4.4). A total of 38 Class Ills are already known 
within the ISOCAM survey area, so that ^ 68 unknown 
Class Ills may remain to be found. Since it was noted in 
Sect. 2.4 that ~ 30 sources detected at 6.7 pm might be 
unidentified cluster members, about half of the missing 
Class Ills may have been actually seen by ISOCAM. 

We also note that a large proportion (80/123) of the 
Class II sources are closely associated with the densest 
part of L1688 (see Fig. 1). Assuming the same propor- 
tion applies to Class Ills, we would expect ^ 44 unknown 
Class III sources to be located within the CS contours of 
Fig. 1. A total of 39 uncl assifi ed ISOCAM sources (also de- 
tected by Barsony et al. 1997) lie within these CS contours 
where the number of detected background stars should 
be small due to high cloud extinction. Most of these 39 
sources might thus be yet unidentified Class III YSOs. 
These candidate Class III sources are listed in order of 
decreasing ajj^^"* in Table 5. 



4. Luminosity estimates and luminosity function 

Most of the new YSOs identified by ISOCAM are weak 
IR sources which were not detected by IRAS and were not 
observed in previous ground-based mid-IR surveys (dedi- 
cated to bright near-IR sources). They likely correspond 
to low-luminosity, low-mass young stars. In Sect. 4.1 be- 
low, we derive stellar luminosity estimates for Class II and 
Class III objects using published near-IR photometry from 
Barsony et al. (1997). In Sect. 4.2, we provide mid-IR es- 
timates of the disk luminosities, Ldisk, for Class II YSOs. 
Finally, calorimetric estimates of the bolometric luminosi- 



ties, £boi, for Class I YSOs are calculated in Sect. 4.3. The 
luminosity function of the p Ophiuchi embedded cluster is 
then assembled in Sect. 4.4. 

4.1. Dereddened J-band stellar luminosities for Class II 
and Class III sources 

The J-band flux provides a good tracer of the stellar lumi- 
nosity for late- type PMS stars (i.e., T Tauri stars) because 
the J-band is close to the maximum of the photospheric 
energy distribution for such cool stars. It is also a good 
compromise between bands too much affected by interstel- 
lar extinction at short wavelengths (very few p Ophiuchi 
YSOs have been detected in the V, R, or I bands), and 
the H, K and mid-IR bands which arc contaminated by 
intrinsic excesses. Greene et al. ( p_994 ) showed that there 
is a good correlation between the dereddened J-band flux 
and the stellar luminosity derived by other methods. They 
pointed out that in p Ophiuchi this correlation is roughly 
consistent with a more theoretically based correlation ex- 
pected for 1- Myr ol d PMS stars following the D'Antona 
& Mazzitelli ( 1994 ) evolutionary tracks. Mo re rec ently, 
Strom et al. ( 1995 ) and Kenyon & Hartmann ( 1995 ) used 
the same model PMS tracks to directly convert dered- 
dened J-band fluxes into stellar masses. We adopt a simi- 
lar approach here. 

4.1.1. Extinction estimates from near-IR colors 

The main difficulty and source of uncertainty with this 
method is due to the foreground extinction affecting the 
J-band fluxes. One must estimate the interstellar extinc- 
tion toward each source in order to correct the observed 
J-band fluxes. We have used the observed near-IR colors 
to estimate the J-band extinction. The (J — H) color ex- 
cess is most suitable for this purpose (e.g. Greene et al. 
1994 ) since the dispersion in the intrinsic (J — iJ) colors of 
CTTSs is s mall a nd observatio nally well determined (cf. 
Strom et al. |1989| , Meyer et al. [19971) . 

The reddening law quoted by Cohen et al. (1981), 



which is determined for the standard CIT system, should 
be ap plicable to the JHK photometry of Barsony et al. 
(|1997|). We have thus used: 



Av = 9.09 X [(J -H)-{J - H)o] 
when the (J — H) color is available, and 
Av = 15.4 X [{H -K)-{H - K)o] 



(1) 



(2) 



otherwise, where we have adopted (J — H)o — 0.85, and 
{ H - K)q = 0.55 from the CTTS results of Meyer et al. 
( |1997D . The dereddened Jo magnitude is then obtained as 
Jo = J - 0.265 X Av, or Ho = H - 0.155 x Av for those 
stars not detected in J (with Av estimated from H-K). 
The dereddened Jo (or Hq) can then be converted into an 
absolute J-band (H-band) magnitude, Mj (Mh), via the 
distance modulus of the p Ophiuchi cloud equal to 5.73 
(for d ^ 140 pc). 
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The uncertainties on Mj and Mh result from the 
typical uncertainties on the J, H, K magnitudes and 
on the intrinsic colors (J — H)q and [H — K)q. With 
cr(J) = 0.1, (t{H) - <t{K) = 0.05 (Barsony et al. |1997D , 
an d a {(J - H)o)= 0.1, (y{{H - K)o)= 0.2 (Meyer et 
al. 1997), we obtain the following typical uncertainties: 
a{Mj) = 0.39 mag, and a{MH) = 0.60 mag. In addition, 
the uncertainty on the cluster distance (140 ± 10 pc) in- 
duces a maximum systematic error of ±0.15 mag on Mj 
and Mh- 

4.1.2. Relationship between Mj and for young 
T Tauri stars 



The absolute J-band magnitude Mj can be directly con- 
verted into a stellar luminosity Li, if the effective stel- 
lar temperature Tcff is known: log]^g(L^) = 1.89 — 0.4 x 
{Mj + BCj{Tcs j), where BCj is the bolometric correc- 
tion for the J band depending only on Teg. Pre- main se- 
quence objects in the mass range 0.1 ^ Mi, ^ 2.5 Mq are 
cool sub-giant stars with typical photospheric tempera- 
tures Toff ~ 2500-550 K ( e.g. Greene & Meyer |l995, 
D'Antona & Mazzitelli |1994 ). In this temperature range 
(0.34 dex wide), the photospheric blackbody peaks close 
to the J band (1.2/im), so that the J-band bolometric 
correction spans only a limited range, 1 ^ BCj ^ 2, cor- 
responding to a total shift in luminosity of only 0.4 dex. 
Therefore, if we use a (geometrical) average value for the 
effective temperature, Toff^ 3700 K, we should not make 
an error larger than ±0.2 dex on log]^Q(L^). 

Furthermore, since bright Afj sources tend to be of 
earlier spectral type than faint sources, we can in fact 
achieve more accurate luminosity estimates. Indeed, PMS 
stars are predicted to lie within a well-defined strip of 
the HR diagram. This is illustrated in Fig. 4 which 
displays model evolution ary t racks and isochrones from 
D'Antona & Mazzitelh ( |1998[ ) on a Mj-logio(rofr) dia- 
gram for PMS stars with ages between 0.1 Myr and 3 Myr . 
We have used the com pilatio ns of Hartigan et al. (1994), 
Kenyon & Hartmann ( |1995D , and Wilking et al. ( |1999D 
to derive an approximate linear interpolation for BCj: 
BCj = 1.65 + 3.0xlogio(3700/roff)- Figure 4 shows that, 
for a given observed value of Mj, the possible range of 
logio(T'off) is reduced to less than ~ 0.15 dex, inducing a 
maximum error of ±23% (±0.09 dex) on Li,. Based on 
Fig. 4, we have adopted a linear relationship between AIj 
and Teff: log(reff/3700) = -0.055 x {Mj - 4.0) (see the 
heavy dashed line in Fig. 4). This leads to the following 
Li, — Mj conversion: 



logio(L,) = 1.49-0.466 x Mj 



(3) 



The uncertainty resulting from this conversion should 
be on the order of 0.045 dex for 0.01 < < 10.0 Lq. We 
therefore estimate that the total uncertainty on logj^gl-^*) 
is fT(logio I-*) = (0.0452-f (0.466 X(t(Mj))2)1/2 = 0.19 dex. 

The Li, — Mj conversion described above cannot be 
applied to sources undetected in the J band. Instead, we 
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Fig. 4. Mj— logio(rcft) "HR" diagram showing the well- 
defined strip of the HR diagram occupied by the PMS 
tracks of D'Antona & Mazzitelh (1998). For Mj = 4.0 for 
instance, we see that logio(Tofr) can span a range of only 
~ 0.15 dex for stellar ages between 0.1 Myr and 3 Myr. We 
thus adopt the Mj— logio(Tcff) linear relationship dis- 
played as a heavy dashed line which corresponds to the 
conversion: logio(L^) = 1.49 - 0.466 x Mj. 

use the H-band magnitude, along with the extinction esti- 
mate derived from the {H — K) color, but with the addi- 
tional complication that the circum stella r (disk) emission 
cannot be neglected. Meyer et al. (1997) found that the 
H-band circumstellar excess is on the order of 20 % of the 
stellar flux, on average, for the Taurus CTTS sample of 
Strom e t al. (1989 ) (for a p Ophiuchi sample, see Greene 
& Lada 1996). This excess, expressed as a veiling index 
rji = Fi^^jFll (e.g. Greene & Meyer |1995D, is equal to 



^ 0.2. Accordingly, we have applied a systematic correc- 
tion (5H = -2.5 logio(l + th) ~ 0.2 mag. The - Mh 
relationship obtained in a way similar to the J-band rela- 
tion is then: 



logio(i*) = 1.26 - 0.477 x {Mh + 0.2) 



(4) 



The total uncertainty on log]^Q(L^) derived from Mh 
is then cr(logioi,) = (0.045^ + (0.477 x cr(Afo))2)i/2 ^ 
0.29 dex. 

This method can also be applied to Class III YSOs, 
using different values for (J — H)q and {H — K)^. We 
have derived Li, for all Class III YSOs using the following 
relationships: Ay = 9.09 x [(J- 77) - 0.6], or Ay = 15.4 x 
[(77-X) -0.15]; and log(L^) = 1.49-0.466 x (J-0.265 x 
Ay - 5.73), or log(L^) = 1.26 - 0.477 x {H - 0.155 x 
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Av - 5.73) (the H-band IR excess for Class III YSOs is 
negligible). 

The resulting Mj, Mh, Ay, and L^, estimates are listed in 
Table 3 for Class II YSOs, and in Table 4 and Table 5 for 
Class III YSOs. 

4.2. Disk luminosities for Class II YSOs 

Since the SED of an embedded Class II YSO peaks in 
the mid-IR range, the ISOCAM fluxes should be approx- 
imately valid tracers of the total, bolomctric luminosities 
(Lboi)- To estimate Lboi for weak Class II sources, an em- 
pirical approach thus consists in using this -F^^^^-^^boi 
relationship after proper calibration on a sub-sample of 
(brighter) objects for which the luminosity can be derived 
by a more direct method. This approach has been adopted 
by, e.g., Olofsson et al. ( 1999 ). Here, we have used the Li, 
estimates of Sect. 4.1 to check that a correlation is actually 
present between L^, and the mid-IR fluxes. Figure 5 dis- 
plays FJ'' '^ (corrected for extinction) as a function of Li, 
for the 104 Class II sources detected both in the near-IR 
and in the mid-IR range. 
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Fig. 5. Correlation between F},'^-'^ (corrected for interstel- 
lar extinction with ^14.3 — 0.03 x Ay) and L*. The ob- 
servations are compared with a purely reprocessing disk 
model with Ldisk = 0.25 x and a stellar contribution 
given by Eq. 6. The continuous line corresponds to an aver- 
age disk inclination [i = 60°), while the dotted and dashed 
lines correspond to more extreme inclinations [i = 0° and 
i = 85°). 

A correlation is found, showing that, despite some scat- 
ter, the ISOCAM fluxes can be used to give rough esti- 
mates of the stellar luminosities of Class II YSOs. This is 
useful for the few ISOCAM sources of our sample which 
have not been detected at near-IR wavelengths. 

The mid-IR emission of Class II YSOs is usually inter- 
preted as arising from warm dust in an optically thick 
circumstellar disk. Using a simplified disk model (e.g. 



Beckwith et al. 1990), it is easy to show that any ob- 
served monochromatic flux in the optically thick, power- 
law range of the disk SED is simply proportional to the 
total disk luminosity Ldisk divided by the projection fac- 
tor cos(z), where i is the disk inclination angle to the 
line of sight. In the Beckwith et al. (1990) model, the 
disk is parameterized by a power-law temperature pro- 
file with three free parameters, Tq, tq, and q, such that: 
r(r) = To X {r/roy. Here, we have adopted Tq = 1500 K, 
meaning that the disk inner radius is at the dust sublima- 
tion temperature, and q — 2/3, corresponding to an IR 
spectral index am = 2/g-4 = -1.0 typical of CTTS 
spectra. 

We must, however, account for the fact that the stellar 
emission itself is not completely negligible in the mid-IR 
bands, especially at 6.7 /im. A simple blackbody emission 
at n = 3700 K gives 

F,*^/ = 0.082 X (L./lLo) Jy (5) 
= 0.021 X (i./lLg) Jy. (6) 

Since -F'Jt^'^ is significantly smaller than F^'l, only Fj'*-^ is 
used to estimate Ldisk- Assuming cos(i) — 0.5 (i = 60°) in 
all cases, we have computed Ldisk for each of the Class II 
YSOs with L^ and Ay estimates from Sect. 4.1, based on 
the following relationship: 

Ldisk/1 Lo = (FJ4-3 X io+0-4x^-3 _ Fj4/)/i.80 Jy, (7) 

where A14.3 = 0.03 x Ay and F^\^ from Eq. 6. The results 
are given in the last column of Table 3. 

According to this model, the mid-IR flux is a direct 
tracer of the disk luminosity, and the F^^^ — L^, cor- 
relation of Fig. 5 simply expresses that Ldisk correlates 
with L^. The origin of Ldisk is either the release of grav- 
itational energy by accretion in the disk, or the absorp- 
tion/reprocessing of stellar photons by the dusty disk. In 
the latter case, Ldisk is naturally proportional to L^,. The 
fraction Ldisk/L* of stellar luminosity reprocessed by the 
disk depends on the spatial distribution of dust. In the 
ideal case of an infin ite, sp atially flat disk, this fraction is 
0.25 (Adams & Shu [19861) . If the disk is flared, Ldisk/L* 
is larger, while it is smaller if the disk has a inner hole. 
The theoretical F^^^ — Li, correlations plotted in Fig. 5 
correspond to Ldisk — 0.25 x Li, and to three represen- 
tative inclination angles. The fact that this simple model 
accounts for the observed correlation quite well, suggests 
that the disks of most p Ophiuchi Class II YSOs are pas- 
sive disks dominated by reprocessing. This is consistent 
with recent estimates of the disk accretion level in Taurus 

1991 ). 



CTTSs (e.g. GuUbring et al. |199g). The typical disk ac- 
cretion rate of a CTTS is estimated to be lO^^Mo/yr, 
corresponding to an accretion luminosity Lace — 0.025 Lq 
for i?^ = 3Ro and Mi, = 0.25 M© (i.e., L^ - 0.25 Lq at 
1 Myr). In this case, the luminosity due to reprocessing is 
^ 5 times larger than the accretion luminosity in the disk. 

On the other hand, 37 sources (among a total of 104) 
are located above the passive disk model lines in Fig. 5. 
These are good candidates for having an active disk with 
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an accretion rate typically larger than ^ 5 x (L*/0.25Lo) x 
lO-SMo/yr. 

Overall, we find that the median Ldisk/i* ratio is 0.41 
for the 93 Class II sources detected in the near-IR and 
with F^'^-^ > ISmJy (i.e. the completeness level derived 
in Sect. 2.4). Using this ratio, we have derived rough esti- 
mates of the stellar luminosities of the 15 Class II YSOs 
which have no near-IR photometry (see Table 3) as fol- 
lows: 

L*(Lq) « idisk/0.41 « (0.97xl.6/0.41)x(^^,i4.3/i 8 jy) (8) 

[The factor 0.97 corresponds to a typical stellar con- 
tribution of 3% at 14.3 /im, while the factor 1.6 is an 
average extinction correction at 14.3 /im (corresponding 
to < Ay >= 17 mag).] 



4.3. Calorimetric luminosities for Class I sources 

The most direct method of estimating the total luminosi- 
ties Lboi of embedded YSOs consists in integrating the 
observed SEDs (cf. WLY89). However, since most of the 
p Ophiuchi Class II and Class III YSOs are deeply embed- 
ded within the cloud {Ay ~ 10), only a negligible fraction 
of their bolometric luminosity can be recovered by finite- 
beam IR observations (e.g. Comeron et al. 1993 ). We thus 
do not attempt to derive calorimetric estimates of Lboi for 
these sources. In contrast, the calorimetric method is be- 
lieved to be appropriate for Class I YSOs since these are 
self-embedded in substantial amounts of circumstellar ma- 
terial which re-radiate locally the absorbed luminosity (cf. 
WLY89 and AM94). Using our new mid-IR measurements, 
we have evaluated the calorimetric luminosities (Lcai) of 
the 16 Class I YSOs observed in our survey. Only 7 of 
them have reliable IRAS fluxes up to 60 or 100 /im (IRS54, 
IRS44, GSS30, IRS43, EL29, IRS48, IRS51). For these, the 
median of the ratio of Lcai(6.7-14.3/im) to Lboi is found 
to be 9.8, suggesting that the typical fraction of a Class I 
source's luminosity radiated between 6.7 and 14.3 /xm is 
~ 10 % . Assuming that this ratio is representative of all 
Class I YSOs, we have derived estimates of Lboi for the re- 
maining 9 weaker Class I sources (i.e., CRBR85, LFAM26, 
LFAMl, WL12, IRS46, CRBR12, IRS67, CRBR42, WL6). 
These luminosities are listed in Table 2. 



4.4. Luminosity function of the p Ophiuchi embedded 
cluster 

Combining the luminosities determined in Sect. 4.1 for 
the sources detected in the near-IR with the estimates 
from Ldisk for the sources without near-IR measurements 
(Sect. 4.2), we have built a luminosity function for Class II 
YSOs which represents a major improvement over previ- 
ous studies (see Fig. 6a). In terms of L^,, the completeness 
level for this population can be estimated from the F^'^-^ 
completeness limit derived in Sect. 2.4 (J^J^ '^ ~ 15mJy) 
using Eq. 8: L^™™P(ClassII) = 0.032 Lq. While the lu- 
minosity function previously published by Greene et al. 



(1994) included only 33 (bright) Class II sources and suf- 
fered from severe incompleteness below L^, ~ 1 — 2 Lq, our 
present completeness level is a factor ^ 30-50 lower. The 
new luminosity function shows a marked flattening in log- 
arithmic units at L^, ^ 2Lq, well above our completeness 
limit. This important new feature is discussed in Sect. 5 
below. 

Based on the Lcai estimates of Sect. 4.3, a new 
bolometric luminosity function for the 16 Class I YSOs 
of p Ophiuchi is displayed in Fig. 6b. The associated 
completeness level is derived from F^'"^ = lOmJy and 
= 15mJy (Sect. 2.4) using Lboi/icai(6.7-14.3/im)« 
9.8 (Sect. 4.3): L^™'°P(ClassI) = 0.017 Lq. The median 
iboi for Class I YSOs is 1.6 L©, which is ~ 8 times larger 
than the median of Class II YSOs (0.20 Lq). The lumi- 
nosities of Class I YSOs span a range of two orders of mag- 
nitude between ~ 0.1 Lq and ~ 10 Lq, which is roughly as 
wide as the luminosity range spanned by Class II YSOs. 
The comparatively large value of iboi for Class I YSOs 
is probably due to a dominant contribution of accretion 
luminosity as expected in the case of protostars. 

In Fig. 6c, we plot the luminosity function of the 55 
Class III sources that are located within the CS con- 
tours of Fig. 1 and for which we have enough near-IR 
data to derive according to the procedure described in 
Sect. 4.1. This sample comprises 19 confirmed Class Ills 
from Table 4 together with 36 candidate Class III sources 
from Table 5. It might be contaminated by a few back- 
ground/foreground sources but is characterized by a rel- 
atively well defined completeness luminosity. From the 
completeness level F^ '^ = lOmJy and using Eq. 5 with 
an average extinction correction corresponding to Ay = 
17 mag, we get L/'""^'(ClassIII) = 0.20 Lq, which is - 6 
times higher than Lj,'^°™''(ClassII). Deep X-ray observa- 
tions with XAIM should improve the completeness lumi- 
nosity for Class III YSOs by an order of magnitude in 
the near future (cf. discussion by Grosso et al. 2000 in a 
companion paper). 



5. Modeling of the luminosity function 

If the mass-luminosity relationship of PMS stars were a 
simple, fixed power-law function, the shape of the lumi- 
nosity function would directly reflect the underlying mass 
function. Unfortunately, the mass-luminosity relation is 
a complex function with inflexion points and is strongly 
age-dependent. Deriving the mass function from the lumi- 
nosity function thus requires knowledge of the stellar age 
distribution. The p Ophiuchi embedded cluster is believed 
to be yo unger than most known star clusters (e.g. Luhman 
& Rieke 1999 ), with a typical age on the or der o f, or less 
than, 1 Myr (e.g. WLY89, Greene & Meyer |l995| ). 

The 123 Class II sources identified with ISOCAM are 
presently the most complete sample of young stars avail- 
able in the cluster. The corresponding mass function is 
estimated down to 0.055 Mq in Sect. 5.1. We investigate 
the effect of including Class III YSOs by modeling the lu- 
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Fig. 6. Luminosity functions (LF) (a) for the 123 Class II YSOs (continuous histogram with statistical error bars). 
The function corresponding to a similar histogram shifted by half the 0.2 dex bin size is shown as a thin curve to 
illustrate the level of statistical fluctuations due to binning. The LF of 33 Class II sources from Greene et al. (1994) is 
displayed as a darker histogram (rebinned to 0.2 dex bins, and rescaled to d = 140 pc for better comparison with the 
new LF). The typical ±1(T uncertainty on the position of each bin resulting from the individual uncertainties on Li, 
is indicated, along with the effect of the distance uncertainty; (b) for the 16 Class I YSOs ; and (c) for the sample of 
55 (19 confirmed + 36 candidate) Class III sources from Table 4 and 5. The LF of (a) is shown in the background for 
comparison. The dashed vertical lines show the respective completeness levels. 



minosity function of the 135 Class II and Class III objects 
located inside the CS contours of Fig. 1 (Sect. 5.3). 

5.1. Constraints on the mass function of Class II YSOs 



We have modeled the observed luminosity function us- 
ing a two-segment power-law mass function with two 
free parameters, Mflat and ai: dA''/dlogM* cx M*"i for 
< Mflat, and diV/dlogM* cx M*-^ '' for > Mflat- 
The high-mass range (for M^, > Afflat) is simply taken to 
be a power-law with index —1.7 (whereas the Salpeter in- 
dex would be —1.35), in agreement with the IMFs favored 
by Kroupa et al. (1993) and Scalo (1998). We assume that 
the cluster forms stars with a fixed IMF, independent of 
time. The evolutionary tracks of D'Antona & Mazzitelli 
( |l998| )n are used to derive a set of mass-luminosity re- 
lations adapted to the p Ophiuchi PMS population, un- 
der the following three hypotheses about the star forma- 
tion history: Constant star formation rate over a) 0.5 Myr 
(Fig. 7a), b) 1 Myr (Fig. 7b), and c) 2 Myr (Fig. 7c). Cases 
b) and c) correspond to the simplest scenarios consistent 
with current observational constraints (e.g. WLY89). Case 
a) mimics a recent "burst" of star formation suggested by 



We here adopt the tracks of D'Antona & Mazzitelli (L998) 



as they provide the most extensive set of PMS models available 
to date, including a tight grid of masses and ages. In the future, 
we plan to check the details of our results using other tracks 
such as those of Baraffe et al. ( |l998| ) or Palla & Stahler ( |l999| ). 



recent studies (e.g. Greene & Meyer 1995 and Luhman & 
Rieke 1999). These three scenarios are roughly represen- 
tative of our current (imperfect) knowledge of the stellar 
age distribution in the cluster. 

The fitting analysis is performed on the 12 (loga- 
rithmic) luminosity bins above the completeness level of 
0.032 Lq and below 10 Lq by varying ai and Mflat. Due 
to the age ambiguity for a particular (M^,, L*) in the 
PMS tracks for M^ > 2 Mq and age < 3 Myr as a re- 
sult of the transition from convective to radiative interi- 
ors (e.g. Fig. 4), our simplified modeling is not valid for 
Li, ^ 10 Lq. The best fit is shown as a heavy curve for 
each of the three assumed age distributions (Fig. 7a, 7b, 
and 7c). The flattening of the luminosity function below 
2 L0 is well reproduced in all three cases. The best-fit 
model in the 1 Myr case (Fig. 7b) is however not quite 
as good as the two others since it predicts a peak at 
Li, ~ 0.4 — 0.5 Lq where the data show a dip. This is re- 
flected by the largest value 1.41 for the reduced xi = -/^r;;; 
(with = 12 fitted data points and n — 2 free parame- 
ters). The model of Fig. 7c reproduces the data somewhat 
better {xi^ = 1.10) as the predicted peak luminosity moves 
down to Li, ^ 0.2 — 0.3 L©. A marginally better fit is ob- 
tained with the 0.5 Myr model (Fig. 7a, xi = 0.66). In this 
case, the model luminosity function has a local maximum 
at L^, ^ 1 Lq as observed and the flattening for Li, ^ 2 Lq 
is particulary well reproduced. 
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Fig. 7. Best model fits (heavy curves) to the luminosity function of Class II sources obtained for a stellar age 
distribution corresponding to a constant formation rate over (a) 0.5 Myr, (b) 1 Myr, and (c) 2 Myr (the first 0.2 Myr 
have been avoided in the PMS tracks to account for a typical protostellar timescale). In (d), similar best fits are 
shown for the combined luminosity function of 135 Class II and Class III sources inside the CS contours of Fig. 1. The 
vertical dashed lines show the completeness levels of the respective luminosity functions (see Fig. 6). 



The best-fit values of Mflat along with the formal 1 a 
errors resulting from the fitting procedure are as follows: 
0.33 ± 0.04 AIq, 0.55 ± 0.09 Mq, and 0.74 ± 0.13 for 
Fig. 7a, 7b, and 7c, respectively. These Mflat values follow 
the expected trend that the older the stars, the higher 
the derived masses for the same luminosities. The final 
error on the determination of Mflat is clearly dominated 
by the uncertainties on the stellar age distribution. We 
adopt the following average value: Mflat = 0.55 ± 0.25 Mq. 
The best-fit values for ai are -0.02 ± 0.15, -0.27 ± 0.13, 
and —0.19 ±0.13, respectively. There is no correlation with 
age in this case. We conservatively adopt ai = —0.15 ± 
0.2. These constraints on the mass function are valid for 
0.032 Lq < Li, < 10 Lq, which approximately corresponds 
to 0.055 Mq < M* < 2M0. 

Under the same three assumptions about the stellar 
age distribution as above, we have also derived a range of 
stellar masses directly from the Mj magnitude for each of 
the 123 Class II YSOs. We display the results as a mass 
function in Fig. 8, where the vertical error bars reflect 
the uncertainties induced by the various age assumptions 
rather than the statistical errors. The best two-component 
power-law mass function derived above is superposed as a 
heavy curve. The error bars are particulary large for the 
two mass bins close to IM©, i.e., just above Mflat. This 
is due to a shift of a large number of stars beyond Mflat 
when the stellar ages are varied upward. The global shape 
of the mass function is however fairly well determined: It 



is basically flat (in logarithmic units) at low masses down 
to 0.055 Mq and shows a steep decline beyond ~ 0.5 Mq. 




M, (Mq) 

Fig. 8. Mass function of the 123 Class II YSOs. The range 

of possible age distributions (see Sect. 5.1) induces a range 
of masses for each star and thus an uncertainty on the de- 
rived histogram displayed as vertical error bars. The best 
two-segment power-law mass function of Sect. 5.1 (heavy 
solid curve), and the effect of binarity as a function of the 
binary fraction from 50 % to 100 % (light curves - see 
Sect. 6.2) are shown. The vertical dashed line marks the 
completeness level of 0.055 Mq. 



14 



S. Bontemps et al.: ISOCAM observations of the p Ophiuchi cloud 



5.2. Possible deuterium feature in the luminosity 
function of Class II YSOs 

Some of the features (peaks and dips) apparent in the 
observed luminosity function might be real as they are 
also predicted by the models (see Fig. 7). The feature ex- 
pected in the luminosity function of young stars as a result 
of deuterium burning during PMS contraction is particu- 
larly interesting since its location is age-dependent (e.g. 
Zinnecker et al. 1993). As already noted, the model fit ob- 



tained under the 0.5 Myr scenario is particularly good as it 
predicts a peak at L^, 1.5 Lq, consistent with the obser- 
vations (see Fig. 7a). In the model, this 1.5 peak corre- 
sponds to an inflexion point in the mass-luminosity rela- 
tion due to deuterium burning in ~ 0.4-Myr old stars with 
~ 0.5 M0. If this peak in the luminosity function is 
confirmed, it will provide strong evidence that the popula- 
tion of Class II YSOs is particularly young in p Ophiuchi. 

5.3. Effect of Class III YSOs 

There is most probably a significant population of 
Class III objects embedded in the cloud (e.g. Sect. 3.5). 
To evaluate how the presence of this population may af- 
fect our conclusions on the IMF of the cluster (see Sect. 6 
below), we here consider the luminosity function of a com- 
bined sample of Class II and Class III sources. The sample 
comprises the 80 Class II and 55 Class III YSOs located 
inside the CS contours of Fig. 1. It should be complete 
down to ^ 0.2 Lq. (Note, however, that only 19 of the 
55 Class III sources of this sample are confirmed YSOs 
at this stage - see Sect. 4.4.) The associated luminosity 
function, displayed in Fig. 7d, is very similar to the lumi- 
nosity function of Class II YSOs. In particular, it exhibits 
a fiattening below ^ 1 — 2Lq. 

The results of a modeling similar to that of Sect. 5.1, 
but performed on only 8 luminosity bins between 0.2 and 
10 Lq, arc shown in Fig. 7d. The best fits found under 
the same three star formation scenarios as in Sect. 5.1 
are good (with xt = 0.64, 0.94, and 1.15, respectively). 
The best-fit values for the parameters Mflat and ai of 
the mass function are as follows: Mflat = 0.32 ± 0.04 Mq, 
0.48 ± 0.08 Mq, 0.66 ± 0.10 Mq, and ai = -0.08 ± 0.45, 
0.00 ± 0.35, -1-0.05 ± 0.19 for star formation durations of 
0.5, 1, and 2 Myr, respectively. Combining these results, 
we obtain Mflat O.5±O.2M0 and ai - 0.0 ±0.4. Within 
the uncertainties these values are essentially identical to 
those found in Sect. 5.1 for the sample of Class II YSOs, 
although they apply to a more limited mass range, ap- 
proximately 0.17Mq < M* < 2M0. 

In contrast to the Class II case, the luminosity function 
of the Class II— Class III combined sample does not show a 
peak at L^, = 1.5 Lq (cf. Fig. 7d). This may be understood 
if the average stellar age of the combined population is 
somewhat larger than that of the population of Class II 
sources (as is expected). 



6. Discussion 

6.1. The initial mass function in p Ophiuchi 

Since the p Ophiuchi molecular cloud is still actively form- 
ing stars, the final IMF of the cluster cannot be directly 
measured. The masses of the stars already formed in the 
cluster provide only a "snapshot" of the local IMF (cf. 
Meyer et al. 200C). However, assuming that the mass dis- 



tribution of formed stars does not change significantly 
with time during the cluster's history (an assumption 
made in the models of Sect. 5.1), any snapshot of the 
mass distribution taken on a large, complete population 
of PMS objects should accurately refiect the end-product 
IMF. 

Compared to previous investigations of the IMF in the 
p Ophiuchi cloud based exclusively on near-IR data (e.g. 
Comeron et al. 1993, Strom et al. 1995, Williams et al. 



1995, Luhman & Rieke 1999), the use of mid-IR photom 



etry has allowed us to consider a much larger sample of 
young stars (see, e.g.. Fig. 6). Since our sub-sample of 
Class II YSOs is fairly large (123 objects) and complete 
down to low luminosities, it provides an excellent oppor- 
tunity to derive improved constraints on the p Ophiuchi 
IMF down to low masses. 

In a statistical sense at least. Class II YSOs are be- 
lieved to represent a specific phase of PMS evolution which 
follows the (Class and Class I) protostellar phases, and 
precedes the Class III phase (see Sect. 1). Due to their 
short lifetime ( ^ 10^ yr), protostars make up only a small 
fraction of a young cluster's population (cf. Fletcher & 
Stabler 1994), and can be neglected in the global mass 
function. Furthermore, in contrast to PMS stars, proto- 
stars have not yet reached their final stellar masses. 

Class III objects are more numerous and should thus 
contribute significantly to the global mass distribution. 
Furthermore, it has been suggested that some YSOs evolve 
quickly to the Class III stage, perhaps as early as the 
"birthline" for PMS stars (e.g. Stabler & Walter 1993), 
and spend virtually no time in the Class II phase (cf. 
Andre et al. 1992, Greene & Meyer 1995). Since such ob- 
jects cannot be identified through IR observations, their 
exact number and mass distribution will not be known 
until the results of deep X-ray (and follow-up) surveys are 
available (see Sect. 3.5). However, providing the (range 
of) evolutionary timescale(s) from Class II to Class HI is 
independent of mass, both classes of PMS objects should 
have identical mass functions. The results of Sect. 5.3 do 
seem to support this view, as they suggest that the mass 
functions of the Class II and Class III samples do not differ 
down to 0.17Mq. 

We therefore conclude (and will assume in the follow- 
ing) that the mass distribution of Class II YSOs deter- 
mined in Sect. 5.1 and shown in Fig. 8 currently represents 
our best estimate of the IMF in the p Ophiuchi embedded 
cluster down to Af^ 0.055 Mq. As discussed in Sect. 6.2 
below, this mass function applies to stellar systems rather 
than single stars. 
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6.2. Influence of binary stars 

A large proportion ( ^ 50%) of field stars are in fact multi- 
ple (e.g. binary) systems (e.g. Duquennoy & Mayor 1991). 
This is also true for young PMS stars, and there is a 
growing body of evidence that the binary fraction is even 
larger in a young PMS cluster like p Ophiuchi than for 
main sequence stars in the field (e.g. Leinert et al. 1993; 
Simon et al. 1995). Since the present study is based on 



ISOCAM/near-IR observations which do not have enough 
angular resolution to separate most of the expected bi- 
naries, a significant population of low-mass stellar com- 
panions are presumably missing from the mass function 
derived above. These low-mass companions are hidden by 
the corresponding primaries. 

To estimate the magnitude of this binary effect, we 
show three simple models in Fig. 8 which assume a pop- 
ulation of hidden secondaries corresponding to a binary 
fraction / of 50%, 75%, and 100 %, respectively. In each 
case, we start from a primary mass function with the two- 
segment power law form derived in Sect. 5.1. We then add 
a population of secondaries with masses distributed uni- 
formly (in logarithmic units) between a minimum mass 
of 0.02 M0 and the mass of the primary. We thus assume 
that the component masses are uncorrelated and drawn 
from essentially the same mass function (cf. Kroupa et al. 
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The global mass functions resulting from addition of 
companions to the primary mass function are shown in 
Fig. 8. It can be seen that the global mass functions are 
similar in form to the original primary mass function. 
Neither the position of the break point (A/fiat) nor the 
slope in the high-mass range are affected by the addition 
of companions. However, the slope in the low-mass range 
(ai) steepens as the binary fraction increases. Indeed, the 
power-law index between 0.055 M0 and A/flat = 0.55 M© 
changes from ai = —0.15 to ai = —0.31 for / = 50%, 

= -0.37 for / = 75%, and ai = -0.42 for / ^ 100%. 

In summary, if we account for uncertainties in the 
binary fraction (50% ^ / ^ 100%), our best estimate of 
the single-star mass function in p Ophiuchi is well de- 
scribed by a two-segment power-law with a low-mass in- 
dex a = —0.35 ± 0.25 down to 0.055 Mq, a high-mass 
index a2 = —1-7, and a break (flattening) occurring at 
Mflat 0.55 ± 0.25 Mq. 

6.3. Luminosity/mass functions in individual cloud 
cores 

Our sample of Class II YSOs is large enough that we can 
study the properties of the four sub-clusters associated 
with the dense cores Oph A, Oph B, Oph EF, and L1689S 
(see Fig. 1). All 123 Class II sources but 5 belong to these 
4 sub-clusters. The luminosity functions of Class II ob- 
jects in the individual sub-clusters are displayed in Fig. 9 
along with the 'best-fit' model of Sect. 5.1 (solid curve). 
They all agree reasonably well in shape with both the to- 
tal luminosity function and the model: all four luminosity 
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Fig. 9. Individual luminosity functions for the four sub- 
clusters associated with the dense cores Oph A, Oph B, 
Oph EF, and L1689S (see Fig. 1). These are displayed 
with 0.4 dex bins and scaled by a factor 1/2 for direct 
comparison with the 0.2 dex bin total luminosity function 
(light grey histogram in the background). The light solid 
curves correspond to the model of Fig. 7a scaled to the 
number of stars in each sub-cluster. 

functions are essentially flat over two orders of magnitude 
in luminosity and appear to have a peak at L^, ^ 1.5 Lq. 
The agreement is particularly good for the sub-cluster 
with the largest number of stars, Oph A (see Fig. 9a), but 
even the smallest sub-cluster, L1689S, tends to reproduce 
the shape of the global luminosity function on a smaller 
scale (Fig. 9d). In contrast, for instance, the luminosity 
function derived for the Chamaeleon I cloud based on 



ISOCAM data (Persi et al. pOOOD differs markedly from the 
p Ophiuchi luminosity functions. It does not exhibit any 
peak at 1.5 Lq and is consistent with an older (~ 3 Myr) 
PMS population having a similar underlying mass func- 
tion (see Kaas & Bontemps 2001). 

The similarity of the individual luminosity functions 
suggests similar distributions of stellar ages and stellar 
masses in each of the four sub-clusters. 

6.4. Global aspects of star formation in p Ophiuchi 

After correction for unresolved binaries (assuming a bi- 
nary fraction / — 75%), the total number of Class II 
sources (including companions) down to 0.055 Mq is 
~ 145. Assuming a Class III to Class II number ratio 
of 19/22 as found by X-ray surveys (Grosso et al. 2000 
- see Sect. 3.5), we infer the presence of ~ 125 Class III 
stars (including associated companions) in the same mass 
range. The typical number ratio of Class Is (plus Class Os) 
to Class lis is 18/123 suggesting an additional ^ 21 
embedded YSOs. Altogether, we therefore estimate that 
there are currently ~ 291 YSOs down to ~ 0.055 M© 
including ^ 19 % of brown dwarfs. The average and me- 
dian masses of these objects are ~ 0.35 M© and ~ 0.20 Mq 
respectively. The total mass of condensed objects (includ- 
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ing brown dwarfs) in the cluster is thus estimated to be 
^^^ciust ^ 291 X 0.35 - 102 Mq. (The brown dwarfs with 
0.055 Mq < < 0.08 M© contribute only - 4 % of this 
mass.) Restricting ourselves to L1688 (thus subtracting 
the ~ 10 % contribution from L1689) whose average radius 
is approximately 0.4 pc (cf. CS contours in Fig. 1), we find 



^L1688 _ 262, 



L1688 



92 Mq, 



L1688 



980 stars/pc^ 



and pLi688 ^ ]yj^ J p^^s ^ -^j^gj-g a,nd are the stellar 
number density and stellar mass (volume) density of the 
cluster, respectively. 

Adopting a conservative value of 2 Myr for the cluster 
age, the total mass of = 102 Mq translates into 

an average star formation rate of 5.1 x 10~^ Mg/yr, corre- 
sponding to one new YSO (of typical mass 0.20 Mq) every 
- 4000 yr. 

Lastly, we can derive the star formation efficiency 
(SFE) in L1688, defined as SFE= M,tar/(Mstar + Mgas). 
The total molecular gas mass, Afgas, of L1688 has been 
estimated to range between 550 M© (from C^^O measure- 
ments - Wilking & Ladafgsl and 1500 M© (from CS(2- 
1) data - Liseau et al. [19951 ). Using mJ'^^^^ = 92 M©, 
we thus get SFE^^'^*^ 6 - 14 %, which is somewhat 
lower than previous estimates (> 22% - WLY89). Note, 
however, that active star formation in L1688 appears to 
be limited to the three sub-clusters/dense cores Oph A , 
Oph EF, and Oph B (see Fig. 1 and Loren et al. |l990|) , 
where the local star formation efficiency is significantly 
higher: SFE ^ 31 %, using a total core mass of 200 M© 



(Loren et al. |1990D 
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Fig. 10. Comparison of the pre-stellar mass spectrum 
measured by MAN98 for 58 protocluster condensations 
(dark histogram with statistical error bars) with the YSO 
mass function derived in Sect. 5.1 for 123 Class II sys- 
tems (light histogram with error bars accounting for both 
statistical and age uncertainties). 



6.5. Comparison with the p Ophiuchi protocluster 
condensations 

In an extensive 1.3 mm dust continuum imaging survey 
of L1688 with the IRAM 30 m telescope (11" resolution), 
MAN98 could identify 58 compact starless condensations. 
Molecular line observations (e.g. Belloche, Andre, & Motte 
^001 ) indicate that the condensations are gravitationally 
bound and thus likely pre-stellar in nature. MAN98 noted 
a remarkable similarity between the mass spectrum of 
these pre-stellar condensations and the IMF of Miller & 



Scalo ( |1979D . 

In Fig. 10, we compare the pre-stellar mass spec- 
trum determined by MAN98 with the mass distribution 
of Class II YSOs derived in Sect. 5.1. (As such, both dis- 
tributions are uncorrected for the presence of close binary 
systems.) It can be seen that there is a good agreement 
in shape between the two mass spectra. A Kolmogorov- 
Smirnov test performed on the corresponding cumulative 
distributions confirms that they are statistically indistin- 
guishable at the 95% confidence level. This supports the 
suggestion of MAN98 that the IMF of embedded clusters 
is primarily determined by cloud fragmentation at the pre- 
stellar stage of star formation. The fact that both the pre- 
stellar and the YSO spectrum of Fig. 10 present a break 
at roughly the same mass 0.5 Mq is quite remarkable. 
A small, global shift of the masses by only ^ 30% upward 
or downward in one of the spectra would make them dif- 
fer at the 2a statistical level. Although in absolute terms, 
both sets of masses are probably uncertain by a factor 
of ^ 2 (due to uncertainties in the 1.3 mm dust opacity 
and in the cluster age, respectively), this suggests that the 
protocluster condensations identified at 1.3 mm may form 
stars/systems with an efficiency larger than ~ 50-70%. 

7. Conclusions 

We have used ISOCAM to survey the p Ophiuchi main 
cloud for embedded YSOs down to a completeness level of 
10 — 15mJy at 6.7 /im and 14.3 /im. Our main findings 
arc as follows: 

1. A total of 425 point sources are detected in the ~ 
0.7 deg^ field covered by the survey, of which 211 are 
seen at both 6.7/Lim and 14.3 /xm. The observed distri- 
bution of the flux density ratio F^'^-^ /F^ "^ is clearly 
bimodal with a gap separating two distinct groups of 
sources: a "red" group corresponding to Class I and 
Class II YSOs with optically thick mid-IR excesses, 
and a "blue" group consisting of Class III YSOs and 
background stars (cf. Fig. 3). 

2. The red group consists of 139 cluster members, 71 of 
them being newly identified YSOs. Based on their mid- 
IR colors from 2 /im to 6.7 /im and 14.3 /tm, essentially 
all of these new YSOs are low-luminosity Class II ob- 
jects. This brings the total number of Class II members 
to 123, a factor of 2 larger than previously known. Only 
50 % of these Class II sources have a near-IR excess 
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large enough to be recognizable in a near-IR color-color 

diagram. 

3. Combining near-IR data from Barsony et al. (1997) 



with our mid-IR photometry, we derive stellar lumi- 
nosities for 123 Class II and 74 Class III YSOs. We also 
estimate bolometric luminosities for 16 Class I objects. 
The corresponding luminosity functions (Fig. 6) are 
complete down to ^ 0.02 Lq, ~ 0.03 Lq, and ~ 0.2 Lq 
for Class I, Class II, and Class III objects, respectively. 
The luminosity function of Class II YSOs is essen- 
tially flat (in logarithmic units) below ~ 2 Lq , with 
a possible peak at 1.5 and dip at 0.5 Lq seen in the 
four sub-regions Oph A, Oph B, Oph EF, and L1689S 
(Fig. 9). 

4. The large proportion of Class II objects observed 
above the completeness level for Class III sources 
(M^ ~ 0.2 Mq) suggests that more than 50% of the 
p Ophiuchi embedded YSO population have optically 
thick circumstellar disks at mid-IR wavelengths. In a 
majority of cases, the luminosities of these disks are 
consistent with pure reprocessing of stellar light. Only 
~ 35% of the Class II objects have excess mid-IR lumi- 
nosities suggestive of substantial disk accretion rates. 

5. The luminosity function of Class II objects is well mod- 
eled by a population of PMS stars with ages in the 
range ~ 0.3—2 Myr and a roughly flat mass distribu- 
tion below Mflat ~ 0.5 Mq, i.e., dA^/dlogM* cx M^-"-^'^ 
down to 0.06 Mq. 

6. Pending the results of a deep X-ray census of Class III 
objects, we argue that the mass distribution of Class II 
YSOs is representative of the emergent IMF of the 
embedded cluster. If we account for the presence of 
unresolved binaries, this emergent mass function is well 
described by a two-segment power law with a low-mass 
index ai — — 0.35±0.25, a high-mass index a2 — —1-7, 
and a transition mass Mflat — 0.55 ± 0.25 Mq (where 
dA''/dlogM* oc Mi,°'). We find no evidence for a sharp 
turnover at low masses down to at least ~ 0.06 Mq. 

7. The shape of the mass function for Class II systems 
is statistically indistinguishable from the mass spec- 
trum determined at 1.3 mm by Motte, Andre, & Neri 
( |l998| ) for the pre-stellar condensations of the proto- 
cluster. This supports the conclusion of these authors 
that the IMF may be primarily determined by frag- 
mentation at the pre-stellar stage of star formation. 
It also suggests that the 1.3 mm protocluster conden- 
sations should form stars/systems with an efficiency 
larger than - 50-70%. 
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Table 1. List of the 212 ISOCAM sources identified as members of the p Ophiuchi cluster 



ISO 


Coordinates(J2000)" 


1/ 


„ b 


rf4.3 


„ b 

^14 3 


Identification*^ 


IR"* 


# 

IT 




5 




\Jy] 




\Jy] 




Class 


1 


16''25™3e!8 


-24°15'51" 


0.157 


0.031 


0.249 


0.023 


IRS2 


II 


2 


16''25'"3g!2 


-24°22'40" 


0.075 


0.014 


0.088 


0.007 


B162538-242238 


nil 


3 


16'' 25^39': 7 


-24°26'37" 


0.197 


0.020 


0.449 


0.013 


IRS3 


nil 


4 


16''25'"4in 


-24° 21' 38" 


0.008 


0.008 


- 


- 


B162541-242138 


III 


5 


16''25'"50!6 


-24° 39' 18" 


0.045 


0.012 


- 


- 


IRSIO 


III 


6 


16''25"'56!3 


-24°20'50" 


0.400 


0.042 


0.662 


0.024 


SR4/IRS12 


II 


7 


16''25'"57!5 


-24°30'34" 


0.049 


0.007 


0.009 


0.006 


GSS20 


III 


8 


16'' 26"" 01= 5 


-24°29'47" 


0.034 


0.010 


- 


- 


B162601-242945 


III 


9 


16''26"'01'!5 


-24°25'20" 


0.025 


0.009 


0.032 


0.007 


SKSl-4 


nil 


10 


16''26'"0S2 


-24°23'40" 


0.507 


0.017 


0.341 


0.024 


DoAr21/GSS23 


III 


11 


16'' 26"" 03! 2 


-24° 17' 49" 


0.025 


0.002 


0.006 


0.003 


VSSG19 


III 


12 


16''26'"04!6 


-24° 17' 52" 


0.005 


0.002 


0.012 


0.003 


B162604-241753 


nil 


13 


16''26'"07!0 


-24°27'28" 


0.112 


0.022 


0.128 


0.015 


B162607-242725 


nil 


14 


16'' 26"" 07! 5 


-24°27'43" 


0.032 


0.007 


- 


- 


B162607-242742 


III 


15 


16'' 26"" 08! 7 


-24° 18' 53" 


0.010 


0.003 


- 


- 


CRBR4 


in? 


16 


16''26'"09!7 


-24°34'12" 


e 


- 


e 


- 


SR3/GSS25 


Ill 


17 


16''26""10!4 


-24°20'58" 


0.386 


0.011 


0.43 


0.017 


GSS26 


II 


18 


16''26'"15!8 


-24°19'23" 


0.026 


0.002 


0.005 


0.003 


SKSl-7 


III? 


19 


16''26""17!0 


-24°22'26" 


0.262 


0.005 


0.414 


0.012 


GSS29 


II 


20 


16''26""17!2 


-24°20'23" 


0.146 


0.003 


0.293 


0.004 


DoAr24/GSS28 


II 


21 




-24°23'49" 


0.246 


0.004 


0.363 


0.004 


CRBR12 


nl 


22 


16''26""18!6 


-24° 17' 11" 


0.010 


0.001 


- 


- 


B162618-241712 


III? 


23 


16''26'"18!8 


-24°26'13" 


0.033 


0.008 


0.027 


0.007 


SKSl-10 


nil 


24 


16'' 26"" 18! 9 


-24°28'22" 


0.673 


0.025 


0.575 


0.024 


VSSGl 


II 


25 


16''26""19!0 


-24°23'07" 


0.051 


0.008 


- 


- 


CRBR17 


III? 


26 


16''26'"19!2 


-24°24'16" 


0.056 


0.006 


0.074 


0.004 


CRBR15 


nil 


27 


16''26'"20!7 


-24°08'48" 


0.030 


0.010 


0.014 


0.006 


WSB28 


III 


28 


16''26""21!0 


-24° 15' 42" 


0.037 


0.004 


0.013 


0.005 


B162621-241544 


III? 


29 


16''26'"2l!5 


-24°23'07" 


4.09 


0.12 


19.1 


0.26 


GSS30/GY6 


I 


30 


16'' 26"" 21! 7 


-24°26'02" 


0.028 


0.007 


0.033 


0.014 


GY5 


nil 


31 


16''26'"21!8 


-24°22'49" 


0.066 


0.013 


0.119 


0.025 


LFAMl 


I 


32 


16''26'"22!2 


-24°44'36" 


0.027 


0.009 


0.028 


0.012 


GY3 


nil 


33 


16''26""22!4 


-24° 24' 08" 


0.009 


0.003 


0.017 


0.003 


GYll 


nil 


34 


16''26'"22!6 


-24°22'56" 


0.052 


0.008 


0.084 


0.016 


GY12 


III 


35 


16''26""23!1 


-24°28'49" 


0.018 


0.008 


0.026 


0.014 


GY15 


nil 


36 


16''26""23!4 


-24°21'02" 


2.59 


0.068 


2.56 


0.066 


GSS31/GY20 


II 


37 


16''26'"23^7 


-24°24'40" 


0.395 


0.005 


0.405 


0.007 


LFAM3/GY21 


II 


38 


16''26'"24!0 


-24°43'09" 


0.174 


0.012 


0.206 


0.017 


DoAr25/GY17 


II 


39 


16''26'"24!1 


-24°24'50" 


2.65 


0.10 


2.16 


0.067 


S2/GY23 


II 


40 


16'' 26"" 24! 2 


-24° 16' 14" 


2.73 


0.027 


3.02 


0.030 


EL24 


II 


41 


16''26""25!3 


-24°24'47" 


0.064 


0.008 


0.049 


0.011 


GY29 


nil 


42 


16''26'"27!7 


-24°26'43" 


0.003 


0.001 


- 


- 


VSSG29/GY37 


III? 


43 


16''26""28!0 


-24°41'51" 


0.029 


0.002 


0.115 


0.006 


GY33 


nil 


44 


16''26""28!5 


-24° 15' 45" 


0.023 


0.004 


0.010 


0.004 


B162628-241543 


III? 


45 


16''26'"29P7 


-24°19'08" 


0.015 


0.002 


0.005 


0.003 


LFAM8/SKS1-19 


Ill 


46 


16'' 26"" 30! 5 


-24°22'59" 


0.347 


0.094 


0.303 


0.065 


VSSG27/GY51 


II 


47 


16'' 26"" 30! 9 


-24°31'07" 


0.012 


0.008 






IRS14/GY54 


III? 


48 


16'' 26"" 34! 2 


-24°23'30" 


e 




e 




S1/GY70^' 


Ill 


49 


16'' 26"" 36! 7 


-24° 18' 10" 


0.007 


0.002 






B162636-241811 


III? 


50 


16''26'"36!8 


-24°19'01" 


0.007 


0.003 






B162636-241902 


III? 


51 


16'' 26"" 36! 9 


-24°15'54" 


0.307 


0.004 


0.636 


0.009 


B162636-241554 


nil 


52 


16''26"'38!0 


-24°23'01" 


0.057 


0.049 


0.066 


0.059 


VSSG4/GY81 


nil 


53 


16''26'"38^6 


-24°23'10" 


0.035 


0.045 


0.046 


0.043 


GY84S 


nil 


54 


16'' 26™ 40! 6 


-24°27'16" 


0.089 


0.003 


0.158 


0.004 


GY91/CRBR42 


I 


55 


16'' 26™ 40! 7 


-24°30'53" 


0.004 


0.002 






IRS16/GY92 
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The absolute ISOCAM coordinates have been determined using the Barsony et al. (1997) database as a reference frame. 
The resulting positional uncertainty is estimated to be ~ 3" (rms). 

The quoted rms uncertainties on F^'"^ and Fl*"'^ comprise the random errors estimated from the temporal and spatial 
noises measured on the data. They do not include the systematic errors due to absolute calibration and aperture correction 
uncertainties. The maximum systematic error is estimated to be less than ~ 15 %. 

For a mo re complete cross-identification of th e p Ophiuchi IR sources, see, e.g.. Table 1 of AM94 and Table 3 of Barsony et 
al. ( 1997 ). "B" refers to Barsony et al. (1997). "ISO" refers to the present survey. Twelve ISO sources are completely new 
sources with respect to published IR surveys; they are named after their ISOCAM coordinates and marked by a *. 
Adopted SED classes for the ISOCAM sources (see text in Sect. 3). "nl" and "nil" indicate newly identified Class I and 
Class II YSOs; a question mark indicates an uncertain classification. 

IS016 = SR3, IS048 = SI, ISO90 = WL22, and IS092 = WL16 are extended sources in the mid-IR range, most probably due 
to PAH-like emission. No attempt has been made to estimate their integrated F^'"^ and Fj* "^ fiuxes. 

No precise mid-IR positio n cou ld be derived for IS048 = S1/GY70. The quoted coordinates correspond to the near-IR 
position of Barsony et al. (1997). 

The position of IS053 is offset by ~ 16" to the south of the near-IR source GY84. This large mid-IR offset is most likely 
due to the presence of extended, structured emission from the bright neighboring source SI. 

IS088 = SR24 is a binary system (SR24S/GY167 and SR24N/GY168; e.g. Greene et al. |l994| ) which is not resolved in the 
ISOCAM images. 
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Table 2. Class I YSOs 



ISO Identification ^"^ Lboi'' 
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1.1 



The ISO number refers to the numbering of Table l(available only in electronic form at http:/ /cdsweb. u-strasbg.fr/ ). 
Lboi is estimated by integrating under the observed SED from 1.2 to 60 fim or 100 /im and extrapolating from 60 pm or 
100 fim to 200 fim with a spectral index of am = -1.0. For the 9 weak sources (ISO # 137, 99, 31, 65, 146, 21, 209, 54 and 
133) without reliable IRAS fluxes, Lboi is estimated using a typical I/boi/i'cai(6.7 — 14.3/im) ratio of 9.8 (see Sect. 4.3). 
Although IS0159 = IRS48 and IS0167 = IRS51 formally lie below our practical Class I-Class II limit (ajR^* = 0.55), they 
are still considered as Class I YSOs here (cf. WLY89 and AM94). 



Table 3. Class II YSOs 



ISO 


Identification 
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3.1 


2.5 


77 


GY152 


0.05 








0.037'= 


0.015 


70 


WL2/GY128 


0.05 




2.3 


38.6 


1.1 


0.27 


165 


GY312 


0.03 


6.6 




14.8 


0.027 


0.064 


85 


CRBR51 


0.03 








0.025'= 


0.010 


175 


GY344 


0.02 


6.5 




17.4 


0.030 


0.074 


26 


CRBR15 


0.01 


6.7 




14.6 


0.022 


0.061 


139 


GY260 


-0.03 




3.7 


40.3 


0.24 


0.058 


37 


LFAM3/GY21 


-0.06 


4.5 




14.5 


0.25 


0.33 


121 


WL20/GY240 


-0.07 


3.3 




16.5 


0.86 


0.67 


51 


B162636-241554 


-0.09 


4.4 




7.6 


0.28 


0.43 


95 


WL1/GY192 


-0.11 


4.9 




20.9 


0.15 


0.16 


122 


IRS36/GY241 


-0.11 




4.9 


34.3 


0.066 


0.025 
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171 


GY323 


-0, 


.12 






4.0 


30.9 


0.19 


0.070 


107 


GY213 


-0, 


.15 


5, 


.6 


- 


19.4 


0.075 


0.089 


76 


GY146 


"0, 


.16 






4.1 


43.5 


0.17 


0.025 


120 


IRS34/GY239 


-0, 


.23 






2.4 


34.8 


1.1 


0.19 


144 


IRS45/GY273 
L1689-IRS5'' 


-0, 


.24 


3, 


.6 


- 


19.5 


0.66 


0.45 


204 


-0, 


.25 


2. 


.6 


- 


12.8 


1.9 


1.4 


17 


GSS26 


-0, 


.30 


3. 


.2 


- 


22.6 


0.95 


0.43 


93 


GY188 


-0, 


.36 


6, 


,3 


- 


20.5 


0.034 


0.020 


98 


GY195 


-0, 


.36 


5, 


,5 


- 


20.6 


0.087 


0.063 


53 


GY84 


-0, 


.39 


5, 


,8 


- 


14.1 


0.060 


0.037 


190 


GY450 


-0. 


.39 


9. 


,2 


- 


8.6 


0.002 


0.008 


23 


SKSl-10 


-0, 


.41 


7. 


,2 


- 


8.2 


0.013 


0.018 


13 


B162607-242725 


-0, 


.41 


4. 


,6 


- 


19.5 


0.22 


0.12 


117 


GY235 


-0, 


.43 


4. 


,9 


- 


9.9 


0.16 


0.12 


212 


L1689-IRS7'' 


-0. 


.44 


2. 


,6 


- 


14.7 


1.8 


0.64 


79 


GY154 


-0. 


.44 






5.5 


24.3 


0.033 


0.013 


24 


VSSGl 


-0, 


.49 


3. 


,2 


- 


17.1 


0.97 


0.49 


3 


IRS3 


-0, 


.50 


4, 


,3 


- 


4.6 


0.31 


0.28 


39 


S2/GY23 


-0, 


.51 


2. 


,0 


- 


11.8 


3.7 


1.6 


140 


GY262 


-0. 


.52 


3. 


,5 


- 


23.7 


0.71 


0.20 


40 


EL24 


-0, 


.54 


1, 


,8 


- 


10.0 


4.5 


2.1 


52 


VSSG4/GY81 


-0, 


.54 


5, 


,2 


- 


17.6 


0.11 


0.058 


154 


GY291 


-0, 


.60 


4, 


,3 


- 


23.0 


0.29 


0.058 


94 


B162703-242007 


-0, 


.61 






6.6 


16.6 


0.010 


0.006 


128 


WL4/GY247 


-0, 


.67 


3. 


,2 


- 


19.5 


1.1 


0.24 


59 


WL7/GY98 


-0, 


.69 


4, 


,1 


- 


27.2 


0.37 


0.061 


84 


WL21/GY164 


-0, 


.70 


7, 


,8 


- 


14.2 


0.007 


0.008 


67 


GSS39/GY116 


-0, 


.72 


3, 


,1 


- 


16.3 


1.1 


0.25 


41 


GY29 


-0, 


.73 


5. 


,1 


- 


19.1 


0.13 


0.043 


88 


SR24N/GY168 


-0, 


.74 


2. 


,5 


- 


7.7 


2.1 


0.76 


35 


GY15 


-0, 


.75 


6, 


,2 


- 


11.5 


0.040 


0.019 


164 


GY310 


-0, 


.75 


6, 


,4 


- 


4.0 


0.032 


0.024 


63 


GY109 


-0, 


.77 


5. 


,8 


- 


14.6 


0.064 


0.022 


151 


GY284 


-0, 


.77 


4. 


,9 


- 


7.6 


0.15 


0.065 


43 


GY33 


-0, 


.78 


4. 


,2 


- 


15.5 


0.36 


0.092 


138 


B162726-241925 


-0, 


,78 


7, 


,9 


- 


11.2 


0.006 


0.004 


36 


GSS31/GY20 


-0. 


.79 


1. 


,5 


- 


6.1 


5.9 


1.6 


177 


GY352 


-0. 


.79 


5. 


,2 


- 


17.7 


0.11 


0.031 


197 


B162821-244246 


-0. 


.81 


6. 


,0 


- 


20.3 


0.050 


0.009 


110 


SR21/VSSG23 


-0, 


.81 


1. 


,9 


- 


3.5 


4.0 


1.7 


166 


GY314 


-0, 


,83 


3, 


,4 


- 


6.4 


0.80 


0.26 


9 


SKSl-4 


-0. 


.85 


6. 


,0 


- 


9.5 


0.051 


0.022 


12 


B162604-241753 


-0. 


,86 


6. 


,6 


- 


13.4 


0.027 


0.009 


155 


GY292 


-0, 


,90 


2, 


,7 


- 


10.8 


1.6 


0.37 


88 


SR24S/GY167 


-0, 


,91 


2, 


,5 


- 


5.9 


2.2 


0.72 


19 


GSS29 


-0, 


,91 


2, 


,9 


- 


9.4 


1.4 


0.28 


115 


WL11/GY229 


-0. 


,92 


6. 


,3 


- 


13.8 


0.037 


0.015 


196 


WSBGO*^ 


-0. 


,92 


4. 


,9 


- 


2.8 


0.16 


0.075 


30 


GY5 


-0, 


,92 


6, 


,3 


- 


2.5 


0.036 


0.019 


176 


GY350 


-0, 


,94 


6, 


,4 


- 


6.4 


0.033 


0.014 


72 


WL18/GY129 


-0, 


,94 


4, 


,7 


- 


10.4 


0.19 


0.061 


163 


IRS49/GY309 


-0. 


,96 


3. 


,1 


- 


10.1 


1.1 


0.23 


193 


B162812-241138 


-0. 


,99 


6. 


,2 




6.2 


0.039 


0.016 


2 


B162538-242238 


-1, 


,00 


4, 


,6 




10.9 


0.22 


0.063 


78 


VSSG5/GY153 


-1, 


,02 


3, 


,8 




19.7 


0.55 


0.056 


86 


IRS26/GY171 


-1. 


,04 


5, 


,2 




19.9 


0.12 


0.015 


105 


WL10/GY211 


-1. 


,05 


3. 


,4 




12.5 


0.78 


0.12 


87 


B162658-241836 


-1. 


,06 


6. 


,0 




14.2 


0.050 


0.011 


160 


B162737- 241756 


-1, 


,08 


7, 


,2 




4.5 


0.014 


0.006 


83 


B162656-241353 


-1. 


,08 


4. 


,2 




11.1 


0.33 


0.058 
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25 



32 


GY3 


-1.09 


6.5 


- 


0.4 


0.029 


0.015 


1 


IRS 2 


-1.09 


3.4 


- 


4.9 


0.79 


0.15 


185 


GY397 


-1.10 


6.1 


- 


4.1 


0.045 


0.016 


6 


SR4/IRS12 


-1.12 


2.8 


- 


1.9 


1.5 


0.37 


142 


VSSG25/GY267 


-1.13 


4.0 


- 


9.8 


0.42 


0.065 


20 


DoAr24/GSS28 


-1.14 


3.5 


- 


1.8 


0.72 


0.16 


178 


GY371 


-1.16 


5.2 


- 


7.1 


0.11 


0.026 


62 


GSS37/GY110 


-1.19 


2.7 


- 


8.5 


1.6 


0.18 


89 


WL14/GY172 


-1.24 


5.8 


- 


17.2 


0.059 


0.006 


194 


B162813-243249 


-1.24 


5.2 


- 


5.9 


0.12 


0.024 


116 


B162713-241818 


-1.27 


3.7 


- 


10.7 


0.57 


0.046 


106 


B162708-241204 


-1.33 


4.5 


- 


8.0 


0.24 


0.027 


102 


GY204 


-1.33 


6.2 


- 


1.5 


0.039 


0.011 


168 


SR9/IRS52 


-1.35 


2.7 


- 


-0.7 


1.6 


0.23 


199 


SR13 


-1.37 


3.4 


- 


-0.9 


0.78 


0.12 


187 


SR10/GY400 


-1.39 


4.1 


- 


-3.2 


0.38 


0.052 


68 


VSS27 


-1.53 


2.6 


- 


5.6 


1.9 


0.13 


172 


GY326 


-1.56 


5.3 


- 


8.5 


0.11 


0.008 


56 


WSB37/GY93 


-1.58 


4.7 


- 


1.2 


0.20 


0.024 


38 


DoAr25/GY17 


-1.58 


3.4 


- 


0.7 


0.83 


0.11 


- 


GY256 


- 


- 


4.5 


34.5 


0.11 


- 


- 


GY257 


- 


- 


4.1 


33.2 


0.17 


- 


90 


WL22/GY174 


- 


- 


- 


- 


29.^ 


- 


92 


WL16/GY182 


- 


- 


- 


- 


44.^ 


- 


207 


IRAS16289-2457 


- 


- 


- 


- 


1.3" 


0.51 


123 


NewS 


- 


- 


- 


- 


0.077'= 


0.031 


150 


NewS 


- 


- 


- 


- 


0.21" 


0.083 


195 


NewS 


- 


- 


- 


- 


0.14^= 


0.058 


200 


Newf 


- 


- 


- 


- 


0.95^= 


0.38 


201 


Newf 










0.025'= 


0.010 


202 


Newf 










0.024'= 


0.010 


203 


NewS 










0.7^ 


0.28 


205 


NewS 










0.17'= 


0.067 


206 


NewS 










0.099'= 


0.040 


208 


Newf 










0.058'= 


0.023 


210 


New» 










0.067'= 


0.027 


211 


New^ 










0.068'= 


0.027 



The sample of 21 Class II sources with —0.05 < a^^^ < 0.55 might contain a significant population of transition objects 
(flat-spectrum objects) between Class I protostars and Class II T Tauri stars. 



The ISO number refers to the numbering in Table 1 (available o nly in electronic form at ^ttp:/ /cdsweb. u-strasbg.fr/ ). 
The J-band flux of IS0127 = GY245 is taken from Greene et al. ( |l994| ). 

IS077= GY152, IS085 — CRBR51, and the last 13 sources are not detected in all near-IR bands, and no reliable Mj or Mh 

can be estimated. L* has been obtained using Eq. 8 (Sect. 4.2). 

ISO204 = L1689-IRS5 and IS0212 = L1689-IRS7 refer to the IR sources listed b y Gre ene et al. ( [l994| ) in L1689. 
IS0196 = WSB60 corresponds to the source B162816-243657 in Barsony et al. ( |l997| ). 

Since ISO90 = WL22 and IS092 = WL16 are young early-type stars (see Sect. 3.1), L* cannot be derived accurately using 
the method described in Sect. 4.1. The quoted luminosities for WL22 and WL16 are taken from WLY89 and Comeron et 
al. (1993), respectively. They have been scaled to d = 140 pc. 



The adopted names and the J2000 coordinates of the completely new IR sources are given in Table 1. 
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Table 4. Class III YSOs 



ISO 


Identification 


2-14 


Mj 


Mh 


Av 










[mag] 


[mag] 


[mag] 


[Lq] 


114 


WL19/GY227 
WL5/GY246'' 
GY12'' 


—0.05 




— 1.0 


73.4 


52. 


125 


-1.02 


- 


-0.7 


59.8 


39. 


34 


-1.06 


3.3 


- 


19.5 


0.89 


58 


WL8/GY96'' 
GY289'' 


-1.11 


1.9 


- 


35.5 


4.2 


152 


-1.33 


3.3 


- 


27.5 


0.90 


198 


SR20 


-1.65 


1.9 


- 


4.0 


4.0 


133 


GY253'' 


-1.73 


3.3 


- 


31.2 


0.85 


45 


LFAM8/SKS1-19'' 
GY156'' 


-1.82 


4.2 


- 


25.2 


0.34 


80 


-2.00 


3.4 


- 


22.5 


0.84 


10 


DoAr21/GSS23'' 


-2.00 


0.7 


- 


6.0 


15. 


27 


WSB28 


-2.06 


4.1 


- 


4.3 


0.36 


149 


B162730-244726 


-2.20 


3.7 


- 


10.3 


0.61 


135 


VSSG22 


-2.23 


3.0 


- 


17.1 


1.2 


64 


YSSGll*" 


-2.24 


3.4 


- 


15.2 


0.83 


73 


VSSG3/GY135 


-2.26 


2.3 


- 


15.7 


2.8 


180 


VSSG14/GY372 


-2.51 


2.2 


- 


5.5 


3.0 


184 


IRS55/GY380 
VSSG19'' 


-2.65 


2.6 


- 


6.0 


1.8 


11 


-2.75 


3.8 


- 


3.9 


0.50 


7 


GSS20 


-2.84 


3.1 


- 


4.7 


1.1 


60 


GYlOl*" 


( 0.00) 


- 


2.0 


55.5 


2.0 


61 


GY103'' 


(-0.45) 


- 


2.6 


48.5 


1.1 


126 


GY248'' 


(-0.82) 


4.2 


- 


25.5 


0.36 


101 


IRS30/GY203'' 


(-0.99) 


- 


2.2 


36.5 


1.6 


14 


B162607- 242742 


(-1.43) 


3.5 


- 


20.6 


0.70 


8 


B162601-242945 


(-1.96) 


3.9 


- 


8.1 


0.47 


183 


GY377'' 


(-1.98) 


4.0 


- 


16.0 


0.40 


157 


GY296'' 


(-2.02) 


5.5 


- 


5.1 


0.081 


4 


B162541-242138 


(-2.05) 


5.8 


- 


6.7 


0.060 


69 


GY122 


(-2.17) 


5.5 




2.7 


0.080 


96 


GY193 


(-2.24) 


4.2 




7.4 


0.34 


97 


GY194 


(-2.35) 


4.1 




9.1 


0.40 


188 


GY410 


(-2.36) 


4.0 




10.2 


0.43 


66 


GY112 


(-2.42) 


4.2 




3.6 


0.35 


130 


SR12/GY250'' 


(-2.48) 


3.4 




1.2 


0.83 


5 


IRS 10 


(-2.78) 


2.7 




5.4 


1.7 




IRS50/GY306 




3.7 




11.5 


0.60 


16 


SR3/GSS25 
S1/GY70'' 










100.'= 


48 










1100.'= 



" The ISO number refers to the numbering in Table 1 (available only in electronic form at ittp://cdsweb. u-strasbg.fr/). 
Class III YSOs located inside the CS contours of Fig. 1 (see S ect. 3 .5 and 4.4). 



For the two B stars IS016 = SR3 and IS048 = S1 (e.g. Elias 
Sect. 4.1. The quoted values are taken from Lada & Wilking ( 




cannot be derived using the method described in 
and have been scaled to d = 140 pc. 
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Table 5. Class III candidates located within the CS contours of Fig. 1 





I(l6iit ificcit ion 


2-14 


IVIj 


IVlH 


4,, 


r 


it 




l^IR ) 


[mag] 


[mag] 


[magj 




191 


GY463 


-1.46 


4.4 


- 


25.2 


0.28 


44 


B162628-241543 


-1.63 


3.9 


- 


21.7 


0.45 


81 


VSSG7/GY157 


-2.13 


2.1 


- 


29.8 


3.4 


91 


VSSG8/GY181 


-2.20 


2.4 


- 


22.9 


2.4 


28 


B162621-241544 


-2.23 


2.9 


- 


15.5 


1.4 


18 


SKSl-7 


-2.38 


3.6 


- 


17.4 


0.64 


25 


CRBR17 


( 0.68) 


8.3 


- 


10.1 


0.004 


153 


GY290 


( 0.41) 


- 


3.7 


41.5 


0.31 


100 


B162705-244013'' 


( 0.12) 


- 


- 


- 


- 


109 


GY215'' 


( 0.05) 


- 


- 


- 


- 


104 


GY207 


(-0.31) 


- 


4.6 


39.5 


0.11 


49 


B162636-241811'' 


(-0.36) 


- 


- 


- 


- 


15 


CRBR4 


(-0.41) 


- 


1.9 


58.9 


2.2 


146 


GY278 


(-0.58) 


- 


2.4 


47.8 


1.3 


136 


GY258 


(-0.72) 


8.2 


- 


10.8 


0.004 


55 


IRS16/GY92 


(-0.80) 


- 


5.0 


29.1 


0.079 


173 


IRS53/GY334 


(-0.81) 


- 


2.0 


40.2 


2.1 


71 


GY130 


(-0.87) 


- 


3.7 


38.8 


0.31 


162 


GY309 


(-1.01) 


- 


3.3 


42.5 


0.47 


50 


B162636-241902 


(-1.06) 


6.6 


- 


14.0 


0.025 


57 


B162641-241801 


(-1.22) 


- 


3.1 


32.3 


0.58 


111 


WL9/GY220 


(-1.27) 


5.5 


- 


21.6 


0.086 


22 


B162618-241712 


(-1.50) 


4.2 


- 


26.8 


0.34 


113 


IRS32/GY228 


(-1.52) 


3.6 


- 


18.5 


0.62 


82 


GY163 


(-1.53) 


- 


4.1 


32.8 


0.19 


131 


GY255 


(-1.57) 


4.2 


- 


21.6 


0.35 


169 


GY322 


(-1.59) 


4.0 


- 


16.3 


0.43 


47 


IRS14/GY54 


(-1.85) 


4.7 


- 


16.6 


0.21 


189 


GY412 


(-1.90) 


5.2 


- 


17.6 


0.12 


181 


GY373 


(-1.91) 


6.4 




7.7 


0.032 


174 


GY346 


(-1.93) 


5.6 




19.5 


0.075 


74 


IRS20/GY143 


(-1.97) 


3.8 




16.0 


0.50 


179 


GY370 


(-1.99) 


5.4 




10.0 


0.094 


148 


GY283 


(-2.09) 


4.3 




13.1 


0.32 


42 


VSSG29/GY37 


(-2.14) 


6.8 




6.5 


0.021 


192 


GY472 


(-2.14) 


5.9 




20.1 


0.054 


186 


GY398 


(-2.20) 


5.2 




5.2 


0.11 


158 


GY297 


(-2.22) 


6.0 




0.7 


0.047 


156 


GY295 


(-2.52) 


4.3 




4.8 


0.32 



" The ISO number refers to the numbering in Table 1 (available only in electronic form at |ittp://cdsweb.u-strasbg.fr/ ). 
^ Av and Li, were not derived for these sources as they were detected only in the K band by Barsony et al. ( 19971). 



